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Abstract 


Italian Version Recenti lavori suggeriscono che i vetri di spin di Heisen¬ 
berg possano appartenere alia stessa classe di universality di vetri strut- 
turali. Trovare nn equivalente su reticolo per questo tipo di modclli 
consentirebbe probabilmente di semplificarne lo studio sia dal punto di 
vista numerico che anahtico, consentendo eventualmente di rispondere 
alle molte domande riguardanti la transizione vetrosa che sono irrisolte 
da decenni. 

I liquidi sottoraffreddati hanno comportamenti peculiari che, nel caso 
in cui 1’analogia tra i due tipi di sistema reggesse, dovremmo riscontrare 
nclla fase paramagnetica dei vetri di spin di Heisenberg. 

E’ con questo obiettivo che nel presente lavoro si affronta in maniera 
sistematica la caratterizzazione della fase paramagnetica. Si cerca di 
comprendere il ruolo del panorama energetico, con uno studio appro- 
fondito delle strutture inerenti (in analogia con i liquidi sottoraffred¬ 
dati dehniamo struttura inerente di una conhgurazione il minimo locale 
dcll’energia piu prossimo ad essa), che richicde lo sviluppo di un nuovo 
algoritmo di ricerca. Indaghiamo anche se si trovano tracce di una tran¬ 
sizione dinamica nella fase paramagnetica. Di fatto sia l’esistenza di uno 
scenario energetico complesso, sia quella di una transizione dinamica, 
sono tratti caratteristici della fisica dei liquidi sottoraffreddati. 



Abstract 


English Version Recent work suggests that Heisenberg spin glasses 
may belong to the same universality class than structural glasses. Indeed, 
finding a lattice equivalent for supercooled liquids would probably allow 
easier numerical and analytical studies, that may help to answer long 
standing questions on the glass transition. Supercooled liquids have many 
peculiar behaviors that should be found in the paramagnetic phase of 
Heisenberg spin glasses if the analogy between the two systems holds. 

It is with this motivation that we undertake a study of the paramag¬ 
netic phase of Heisenberg spin glasses. We shall emphasize the role of the 
energy landscape, with a detailed study of the properties of the inherent 
structures (by analogy with supercooled liquids, we name inherent struc¬ 
ture the local minimum of the energy function which is closest to the 
current spin configuration). Finding inherent structures will require the 
development of a new search algorithm. We shall investigate as well the 
existence of a dynamic ’’phase transition” in the paramagnetic phase. As 
a matter of fact, both the existence of a complex energy landscape, as 
well as the existence of a dynamic transition, are distinguished features 
of the Physics of supercooled liquids. 
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Introduction 


Spin and structural glasses are two completely different types of sys¬ 
tems. Nevertheless, deep analogies can be found if one considers some 
dynamic properties, metastability, and the role of competing interactions 
in both types of systems. The main aim of this thesis is to try to give a 
good starting point for a solid correspondence between those two types 
of models. A similar job has been done with non-disordered systems, by 
making a one-to-one mapping between the Ising model and lattice gases. 
Although they represent systems that apparently have nothing in com¬ 
mon, those two models have the exact same behavior if we recognize the 
correct relation between the quantities that appear in each model. 

Establishing a similar analogy between spin and structural glasses 
would probably be mutually beneficial for their study. In particular for 
spin models it is easier to demonstrate theorems, that could extend also 
to structural glasses and supercooled liquids. Besides, numerous approx¬ 
imation methods have been developed for spin glasses, and the numerical 
study of spin glasses is far more manageable than the one of structural 
glasses, since it is possible to apply very performing Monte Carlo algo¬ 
rithms on lattices, so that finding an eventual symmetry relating the 
former with the latter could make us extend the results of one model to 
the other, and vice versa. 

Although it is several decades that supercooled liquids and structural 
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Introduction 


glasses are studied [1], their comprehension is far from being complete. 
In particular it is still not known if there exists a real thermodynamic 
glassy phase at low temperature. A salient feature is that lowering the 
temperature there is a very steep growth of the characteristic time scales 
t of the system. The experimental approach is to define T g as the temper¬ 
ature at which r ~ 10 2 — 10 3 s, but this temperature has only a practical 
meaning. From the theoretical side, the Mode Coupling Theory (MCT) 
[2] offers a framework to describe quantitatively this dynamic slowdown. 
In fact, Mode Coupling Theory predicts a dynamic transition where the 
self-correlation times diverge with power law. Nowadays, the Mode Cou¬ 
pling Transition is known to be a crossover rather than a real transition. 
Nevertheless, the functional forms predicted by MCT are routinely used 
to fit experimental results when the autocorrelation times are in the range 
10 -13 s < r < 10 -3 s. 

Recently, a topological interpretation has been given for the dynamic 
transition predicted by the MCT in which the role of the energy land¬ 
scape, i.e. the topological aspect of the potential energy, assumes a de¬ 
terminant role in the dynamics of these systems. At high temperatures 
the typical configurations in which the system finds itself are closer to 
saddle points of the potential energy than minima, so there is always a 
negative direction in its hessian that permits to the system a fast move¬ 
ment in the configuration space. Decreasing the temperature the typical 
number of negative directions of these saddles starts falling too, so that 
the relaxation times increases as the system has always more difficulties 
in finding these more convenient directions. There is a temperature in 
which the density of negative directions of the hessian goes to zero. This 
change in the topology of the energy landscape is called topological tran¬ 
sition, and it causes the dynamic transition, i.e. a high slowing down 
of the dynamics (but with no divergence and no real phase transition). 
Mode Coupling Theory predicts that the relaxation times diverge with 
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power law at this temperature. This is not what happens because when 
the saddles disappear we enter another regime of dynamics in which the 
system has to ‘jump’ from one minimum to the other, in which the typ¬ 
ical times are much longer since they scale with the exponential of the 
energy barrier. We can so appreciate the divergence predicted by Mode 
Coupling Theory only until its characteristic time scales become com¬ 
parable with those of activated dynamics. Another interesting feature 
that follows from the energy landscape scenario is the concept of inher¬ 
ent structure. It is possible to associate, to any state, the configuration 
corresponding to the local minimum where the system would find itself if 
the energy were to be lowered abruptly in some way. We call that config¬ 
uration inherent structure. What is seen is that at high temperatures the 
energies of the inherent structures almost does not depend on T. On the 
contrary, at lower temperatures these energy suddenly start to decrease 
fairly quickly. The interpretation of this phenomenon is that there is an 
exponentially high number of minima at a certain level, and few with 
lower energy, so as long as we are at high temperature the probability of 
finding an exponentially numerous minimum is almost one, and to feel 
the presence of the lower energy inherent structures we need to be at 
lower temperatures. 

On the other hand, we have a stronger hold on the physics of spin- 
glasses (the second partner in our intended analogy), specially for Ising 
spin glasses. Spin glasses [3, 4] are disordered magnetic alloys, such as 
Mii^Cui-a,, where magnetic moments (localized on the Mn atoms, typi¬ 
cally) interact. If these magnetic moments may point basically anywhere, 
their interaction depending only on their relative orientation, one speaks 
of Heisenberg spin glasses. Yet, if the lattice anisotropies enforce the 
magnetic moments to point mostly to the North or South poles, we have 
Ising spin glasses. 

It is since the beginning of the 90’s that there is clear evidence of a 
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thermodynamic spin glass phase at low temperatures for systems with 
Ising spins and finite-range interactions. Also, there is no dynamic tran¬ 
sition, in the sense that there isn’t any temperature in the paramagnetic 
phase in which the relaxation times grow considerably. Moreover, Parisi’s 
solution of the mean field model yields an infinite-step replica symmetry 
breaking, that implies a continuous sequence of phase transitions once 
the threshold of Tsg is passed, with a continuous fragmentation of the 
valleys of the free energy, and a macroscopically large number of differ¬ 
ent configurations that characterize the glassy phase. The slow down of 
the dynamics in the glassy phase would be due to this complexity of the 
phase space. It appears clear that Ising spins don’t offer the possibility 
of a correspondence with structural glasses. 

In spite of the fact that they are a much more reasonable representa¬ 
tion of real spin glasses, extremely little is known about Heisenberg spins 
glasses. Although it is assumed that a phase transition takes place at 
finite temperature, its nature is not clear yet, since there is a degree of 
freedom that was absent in Ising spins, called chirality, that represents 
the handedness of the non-collinear orientation, and that could yield an¬ 
other phase transition. It is reason of intense debate whether spin and 
chirality transitions are two distinct phenomena, or if they take place 
at the same temperature. Not much more has been investigated, so it 
is worth to make an investigation to see if Heisenberg spin glasses have 
features in common with structural glasses. 

There are some reasons to conceive an analogy between structural 
glasses and some types of spin-glasses. In 1987, Kirkpatrick and Thri- 
umalai [5] pointed out the formal identity between the Mode Coupling 
equations and the exact equations for the dynamics of a particular type 
of spin-glasses, with p-body interactions (the so called p-spin models). 
Unfortunately, this analogy has been established so far only for models 
with long range interactions, where the mean field approximation is ac- 


4 



Introduction 


tually exact. Inspired by this intuition, in the year 2000 G.Parisi and 
M.Mezard [6] made a replica study of structural glasses for finite-size 
systems that suggested that the structural glass transition belongs to 
the same Universality Class of the systems with one step of replica sym¬ 
metry breaking. That means that the order parameter’s distribution is 
not continuous as in the Ising spin glasses, but it is takes only two val¬ 
ues L , because there is only one step of fragmentation of the free energy 
landscape. Unfortunately, a simple model, with short range interactions, 
undergoing a 1-RSB glass transition, is yet to be found. 

There are a couple of clues that make us think that finite-range in¬ 
teraction Heisenberg spin glasses are a good candidate for the corre¬ 
spondence with spin glasses that we have exposed. The first is that 
H.Kawamura in 2005 [7] found a double-peak distribution of the order 
parameter typical of the one-step replica symmetry breaking, so the uni¬ 
versality class may be the same. The second is that up to our knowledge 
the only dynamic study of the paramagnetic phase done up to now [8] es¬ 
timates the phase transition at a temperature 50% higher than its actual 
position, and in an optimistic scenario this overestimation could be ex¬ 
plained if what they found was actually the dynamic transition predicted 
by Mode Coupling Theory. 

The above mentioned clues clearly do not suffice to establish our 
seeked analogy. Hence, if one wants to have stronger evidence of the 
similarities between the two models it is necessary to make a systematic 
research of the aspects under which we would like to remark the analogy. 
We conducted our research under two main aspects: 

1. Since the Mode Coupling Transition is best understood through the 
topological properties of the energy landscape, it is absolutely nat¬ 
ural to start our investigation concentrating ourselves on the energy 

Hn absence of spacial symmetries we have two peaks: one at zero, and one at a 
finite value qea- If, the system is symmetric under reflection there are three peaks 
because we have the reflection of QeAj —qEA 
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landscape of the Heisenberg spin glasses, so we disposed ourselves 
to a systematic research of inherent structures. This type of in¬ 
vestigation has never been done, because Ising spins are discrete 
variables, and Heisenberg spin glasses are still a largely unexplored 
topic. The only energy-lowering routine used in bibliography is the 
quench, which is known to be very badly performing. Hence we had 
to find a better performing algorithm to find inherent structures, 
and we made a deep study of their spatial correlation functions. We 
succeeded in hireling a fast algorithm that permitted us to make our 
research with very large samples, compatible with the thermody¬ 
namic limit, and we noticed that the inherent structures are indeed 
relevant for the physics. First, we found the anharmonic term was 
rather small, so we can explain the energy of a thermalized state 
as harmonic fluctuations around a minimum. Second, lowering the 
temperature the correlation lengths of thermal states and the rela¬ 
tive inherent structures converge. 

2. Mode Coupling Theory predicts a divergence of the relaxation times 
in the paramagnetic phase, so we investigated the equilibrium dy¬ 
namics of the Heisenberg spin glass model to see if we would en¬ 
counter the same behavior. In conformity with the former point, 
we wanted to be in the thermodynamic limit, so we had to use 
non-physical algorithms to obtain thermalized configurations, and 
then we switched to physical dynamics. We wanted to verify if the 
relaxation time admitted a power law scaling for 3-4 decades that 
predicts a transition significantly over the critical temperature, as 
it happens in structural glasses. Also this analysis was successful: 
we found a scaling with exponents of the same order, in the same 
range of times, and the estimate of the possible dynamic transition 
was consistently over the spin glass transition, and two standard 
deviations under the previous estimate given in [8]. 
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After these results the analogy between the two systems looks defi¬ 
nitely more solid. Stronger data for this hypothesis could be given by 
an analysis of the variation with temperature of the density of negative 
directions in the hessian, which is a work actually in progress, and of the 
configurational entropy. 
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1.1 Defining a spin glass model 


1.1 Defining a spin glass model 

In the last decades particular attention has been given to a set of 
models which were meant to describe random magnets with co-presence 
and competition of ferromagnetic and anti-ferromagnetic interactions. 

From the pure modclisation point of view a spin glass has two key 
ingredients: disorder and frustration. A simple example is given in fig¬ 
ure 1.1, where we have 4 Ising spins (they can point only up or down) 
interacting in a loop. The interactions between the spins are randomly 
ferromagnetic or anti-ferromagnetic, so it can happen, as in this case, 
that it is impossible to satisfy all constraints. It is clear that in a very 
large system finding the configuration of spins that satisfies the most 
number of constraints is a very hard task. 



1 


Figure 1.1: A toy model which shows frustration. If the interaction on the bond is 
a the spins want to be parallel and if it is a they want to be 

anti-parallel. The model is frustrated because it is impossible to satisfy 
all the conditions simultaneously. In a spin lattice any loop of spins, in 
which the product of the interactions is negative, is frustrated. 


The first spin glass model was presented in 1975 by S.F. Edwards and 
P.W. Anderson [9] to modelize a dilute solution of Mn in Cu, in which the 
interactions between particles can be ferromagnetic or anti-ferromagnetic 
depending on their mutual distance. 

In the Edwards-Anderson (EA) model the spins are placed in fixed 
positions in a D -dimensional regular lattice. This is because the typical 
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1.1 Defining a spin glass model 


time scales of the movement of the particles is overwhelmingly greater, 
and non-trivial interesting effects can be seen on a regular lattice. 

The model’s Hamiltonian is 

H ^ ^ JxySx ' Syi ( 1 - 1 ) 

<xy> 

where the brackets < • > indicate that the sum is performed only over 
the nearest neighbours. The spins sg are m-dimensional vectors with 
||sj|| = 1. The type of spins has a different name depending on m 's 
value: 

• m — 1: Ising spins 

• m = 2: XY spins 

• rn — 3: Heisenberg spins. 

It is clear that the choice of m depends strictly on the anisotropy of the 
material we want to study. The couplings J xy are random variables with 
zero mean to avoid any bias towards ferromagnetism or antiferromag¬ 
netism, and their standard deviation has to be finite. If we call [■] the 
average over the different realizations over the disorder these conditions 
express as 

[v»] = 0, I4] 1/2 = /(=!)■ (1.2) 

The most used probability distribution functions for the J xy s are the 
bimodal and the Gaussian distribution. We stress that in these systems 
disorder is not only due, like in most models, to thermal fluctuations, 
but it also appears explicitly in the Hamiltonian, giving two levels of 
disorder. We decided to write the averages over the different realizations 
of the couplings with square brackets [•], and the thermal averages will 
be written in angular brackets < • >. 

Experimentally we can find spin glasses in many different systems: 
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• Metals: Magnetic alloys in which we have atoms with strong mutual 
interactions diluted in a non-magnetic metal such as Mn in Cu have 
two-site interactions 


Jxy ^ 


cos(2A ;|\x — y \|) 


\x 


(1.3) 


which depend on the distance and evidently yield the possibility of 
having both ferromagnetic and anti-ferromagnetic interactions. 


• Insulators: An example is Feo.sMno.sTiOa, which comprises hexag¬ 
onal layers. The spins align perpendicular to the layers. Within a 
layer the spins in pure FeTi 03 are ferromagnetically coupled while 
spins in MnTiOs are antiferromagnetically coupled. Hence the mix¬ 
ture gives an Ising spin glass with short range interactions. 


• Non-conventional spin glasses: Spin glass theory has been used 
to study a numerous amount of problems outside the condensed 
matter domain, such as: 


— Optimization problems in computer science 

— Protein folding 

— Polymer glasses, foams... 


1.2 Order Parameters 


A phase transition presents often a spontaneous symmetry breaking 
in the system. The order parameter is the observable to whom we impute 
the breaking of the symmetry. It usually transforms, under the symmetry 
changes, in a homogeneous non-trivial way, so that, when it is zero, 
symmetry is preserved and, when it is non-zero, it is broken. 

It is now assumed, after consistent experimental, theoretical and nu¬ 
merical work, that EA spin glasses have two thermodynamic phases: the 
paramagnetic phase and the spin glass phase. This issue has been quite 
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controversial for years, because the static aspect of a configuration in the 
spin glass phase was not easily distinguishable from a paramagnetic one. 
As a matter of fact the order parameter of a spin glass is not a trivial 
object, and it deserves particular attention. 

1.2.1 Ising Spins 

To make the order parameter definition easier we rely on the compar¬ 
ison with the order parameters in the non-disordered counterparts of our 
models. Thus, for Ising spins we will show how the order parameter is 
defined in the Ising model, and see how to extend this definition to the 
EA Ising model. 

Symmetries and Order Parameter in the Ising Model 

The Ising model’s Hamiltonian is 

H = - s s s$, (1.4) 

<xy> 

with sg — ±1. 

It is immediate to see that its symmetry group is Z- 2 , since making 
the symmetry transformation 

Z-2 '■ {•%} —> {—s#} (1-5) 

the Hamiltonian and the functional integration measure are invariant. 
The order parameter is the magnetization 

X 

which is not invariant under our Z 2 symmetry 

Z 2 : m —> —m (1.7) 

so in the disordered phase, for T > T c , we have m = 0, and for T < T c , in 
the ordered phase, we have m ^ 0. Unfortunately the symmetry break¬ 
ing, hence the phase transition, occurs only in the thermodynamic limit 
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L —» oo because in finite systems the average magnetization is zero also 
in the ordered phase, since the energy barriers that cause the ergodicity 
breaking are infinite only in an infinitely large system. Therefore, if we 
want to extract infinite-size properties from a finite-size system, as it is 
done in numerical simulations, we should better study invariant quanti¬ 
ties that are not negatively influenced by the absence of phase transition. 

The correlation function 


G(x,y) =< s s s$ > 

is clearly invariant under our symmetry transformation 

^2 • y ( $x ) ( ^y) ^X^y 


( 1 . 8 ) 


(1.9) 


The correlation function is strictly binded to the order parameter, 
since 

<m > = < (yy!^)(^y! g i() > — M ^2 X/ < sgs y > ■ 


x y x y 

If the system has periodic boundary conditions the correlation function 
is also translationally invariant so G(x,y) = G(x — y), hence 


1 


< m ‘>= xT, <s ° s <’ > = jv 


1 E G < r T = ir(£= 5 ) t 1 - 11 ) 


where the hat indicates the function’s Fourier transform. We see then 
that the correlation function contains a lot of interesting information over 
the thermodynamic system. 

Symmetries and Order Parameter in the EA Ising Model 

Introducing disorder in the Ising model we have 

H = - y JxySgSg (1-12) 


<xy> 


which, as we said in section 1.1, puts us in the position to have to define 
two levels of disorder. We have two types of symmetry now. One is 
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the global Z 2 related to the positions of the spins belonging to a single 
sample of which we just talked. The second is a local symmetry (gauge 
symmetry) related to the realizations of the disorder. Speaking more 
clearly if we associate to each site x a variable rjg = ±1 it is clear that 
the transformation 


Jxy J X y — VxJxyVy 

S x S x = Vx s x 

keeps invariant the relevant magnitudes H and Z . Notice that the prob¬ 
ability distribution function of the J’s does not change with the gauge 
transformation. 

While the global invariance is verified only in the thermodynamic 
limit, the local symmetry is always true. Let’s see how it is possible to 
study a finite size system. 

Since we want to study the properties of the system no matter the 
disorder realization, we can no longer use s x s y as observable for its in¬ 
vestigation because it is not gauge-invariant. 

To construct a gauge-invariant observable that preserves all the sym¬ 
metries that we need, we resort to a simple trick, redefining the Hamil¬ 
tonian in a harmless way, in order for it to represent two replicas (a) and 
(6) of the same realization of the disorder: 

H m = «, •4,(4“ > 4 o) + 4‘ ) 4 t) ) (1-16) 

<x,y> 

The local symmetry remains the same, while the global symmetry group 
changes then to Z 2 ® Z 2 because the invariance holds also if we flip all 
the spins of only one of the two replicas. 

We can now construct a gauge-invariant object 

q., = 4“ ) 4‘ ) (1.17) 


(1.13) 

(1.14) 

(1.15) 
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1.2 Order Parameters 


that is not invariant under the global symmetry. If we define the Edwards- 
Anderson parameter as 

« = ( L18 ) 

X 

we see that it is a good order parameter since under the global symmetry 

q —* ±q. 

As we did before we relate the square order parameter with the cor¬ 
relation function, which in presence of quenched disorder and translation 
invariance is 

G j{x^y^ Jxy) [<~ QxQy ->] G j{x |/, Jxy) i (1-19) 

SO 

[<^ 2> ] = = = w x - (L20) 

x y r 

1.2.2 Heisenberg Spins 

We have seen how to define the order parameter in the EA Ising model 
and how it is related to some other useful observable. Is there a great 
difference when we pass from m — 1 spin dimensions to m — 3 ? 

Indeed we know that the universality class will not be the same, so 
critical exponents and scaling functions will not be the same. Neverthe¬ 
less those are not the only important differences between the two models. 
It will be in the next chapter that we will focus on the recent results that 
stress the diversity implied by m = 3. 

Here we will concentrate on the symmetries implied by 3-dimensional 
spins, showing how it is possible to identify a new order parameter that 
possibly changes the whole physics of the EA model. 

Symmetries and Order parameter in the Heisenberg Model 

In the Heisenberg model the Hamiltonian is slightly modified 

H = - s x -s y (1.21) 

<x,y> 
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where now the spins s x are 3-dimensional vectors with the constraint 
s 1 = 1. A Heisenberg system is invariant under the rotation or inversion 
of all the spins, so its global symmetry group is 

0(3): s x ^Ts x , TT l = I (1.22) 

and the order parameter 

(1.23) 

X 

is vectorial. As in the Ising case, for a finite system < m >= 0 at all 
temperatures, so there is no point in computing it, and we have to find 
a non trivial quantity in finite volumes, as the correlation function 

G(x, y) =< s x - s y > (1.24) 

that is related to the order parameter in the usual way 

<r n 2> =^J2 G ^ ( L25 ) 

r 

Symmetries and Order parameter in the EA Heisenberg Model 

Knowing already the consequences of the presence of disorder we ap¬ 
ply directly the trick seen in section 1.2.1 and we write the two-replica 
Hamiltonian 

H = -J2 ‘W4“ ) ■ 4“’ + 4" 1 '4 4) ) (1-26) 

<x,y> 

that makes us possible to get the Edwards-Anderson parameter. The 
symmetries are 

• Disorder: Z 2 (local) 

• Thermal: 0(3) <g) 0(3) (global) 
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Our gauge invariant this time is a tensor that we call spin-glass site 
overlap 


r <o(4 = 4t«4,{n 




where a,/3 = 1,2,3 are the components of the spin vectors, thus the 
Edwards-Anderson parameter in the EA Heisenberg model is 


Qa/3 ^ ^ T a/?(^), 

X 

that is related with the correlation function 


G(x,y) = [< tr(T(x)rU : y)) > 


(1 28) 


(1 29) 


in the usual way 

[<«r(QQt)>] = 0) = j, X (1.30) 

r r 

Chirality The order parameter scenario in the EA Heisenberg model 
even more interesting than what we have just exposed. In fact one can 
admit the possibility of breaking the global symmetry Z 2 keeping intact 
the 50(3) symmetry [10]. 

We define chirality of a site along the direction n — 1, 2, 3 as 


ft } 1 (X) Sg—e^ ‘ ($x X Sg- |-e M ) (1-31) 

where is the unit vector in the direction fi. If we make a mirror 
exchange of all the particles respect to a plane the chirality of the sites 
remains the same except for the spins on the plane. A ferromagnetic 
system would not feel the reflection because the system is homogeneous 
and non zero chirality is generated only by thermal fluctuations, and 
not by the energy landscape, so it would not be a good observable. On 
the other side spin glasses are not homogeneous, and the in the ground 
states the spins are not aligned, so they are not invariant under this 
transformation. 
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In compliance with what we did above, we define our chiral gauge 
invariant 

e»(x) = k^\x)k^\x), (1.32) 

from which descends the chiral-glass (CG) overlap in the direction /i, 
which is our order parameter, 

( L33 > 

X 

and we remark its connection with the simmetrized correlation function 
ncc,^ _ [< <?{*)?{$) > + < >] 

\ x y> — 2 

through the usual relation 

[< KE„ >] = ^ £ C™(r) = ±xco (1-35) 

r 

Overlaps A very handy feature of disordered systems is the overlap, a 
magnitude that tells us how similar two configurations are, being equal 
to zero if there is no resemblance, and far from zero in the opposite 
case. There are numerous ways to define an overlap, but all of them 
are not invariant under the global symmetries of a Heisenberg spin glass. 
This could be a problem in numerical simulations, because we could 
think that two identical replicas, who differ only by a rotation, have very 
different configurations. Solutions to this would be very computationally 
expensive, so a symmetry-invariant definition of the overlap would be 
very useful. 

For the chirality the solution is direct, since it is a scalar object and it 
is always invariant, so we can use the chiral glass overlap introduced 
just above. 

For the spin-glass (SG) overlap, for example we could use the Edwards- 
Anderson parameter, but it is not 0(3) invariant. If we look at the 
previous paragraphs we can remark that the square overlap 

Q 2 = tr{r(x)T\x)) (1.36) 
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not only is an invariant, but it also gives an idea of the similarity be¬ 
tween two configurations. It is zero for completely different replicas, and 
about 1/3 when we compute it for a replica with itself and its spins are 
isotropically distributed (in presence of anisotropy it tends to 1). 
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2.1 Ising vs Heisenberg 


2.1 Ising vs Heisenberg 

The original version of the EA model was proposed for Heisenberg 
spins, as it is expected that they represent a good approximation for most 
of the spin glass material, in particular those who posses weak magnetic 
anisotropy. In fact Heisenberg spins are much more present in nature 
than Ising spins, so it is logical to be more interested in them.Nevertheless 
Ising spins have been studied far more because of their greater simplicity. 
It has been argued, in the initial studies of the EA Ising model, that their 
weak magnetic anisotropy present in real magnets causes a crossover from 
the isotropic behavior to the anisotropic, shifting the universality class of 
the systems, so that Ising spins, corresponding to the strong anisotropy 
limit would describe also the behavior of magnets with weak anisotropy 
[ 11 ]. 

2.1.1 Ising ordering in zero field 

It hasn’t been clear whether the three dimensional EA Ising model 
presented a finite-temperature spin-glass transition until in 1985 Monte 
Carlo simulations from Ogielski, Bhatt and Young [12, 13] found solid 
evidence of a finite Tsg■ Later investigations on the critical exponents 
found values that were difficult to compare with an experimental realiza¬ 
tion of the Ising spin glass because at present we know of only one type of 
Ising spin glass really existing in nature that we can study in laboratory, 
an insulating spin glass magnet of Feo.sMno.sTiOs that we have already 
mentioned in section 1.1. We stress that when estimates of the critical 
exponents started to be found the argument mentioned above, that Ising 
spins may fairly represent a weakly anisotropic Heisenberg spin glass, 
collapsed immediately, since they were very far from being compatible. 

A point that raises a lot of controversies is the nature of the SG 
ordered state. There are two dominant visions: the “droplet picture” 
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that claims that EA spin glasses are disguised ferromagnets without a 
spontaneous Replica Symmetry Breaking (RSB) and only two ordered 
states [14], and the “hierarchical RSB picture”, inspired on Parisi’s exact 
solution of the mean held model, which claims that the ordered state 
is accompanied with a hierarchical RSB [15]. In figure 2.1 we show the 
distribution of the order parameter in one case and in the other. 




Figure 2.1: The overlap distribution P(q) in the thermodynamic limit for a finite- 
range spin glass. In figure (a) we show the droplet picture: the spin glass 
is a disguised ferromagnet and has only two ordered states, represented 
with two delta functions at ±qEA- In figure ( b ) we see how with hier¬ 
archical RSB P(q) possesses an additional plateau connecting the two 
delta-function peaks. 

2.1.2 Heisenberg ordering in zero field and chirality 
scenario 

Early numerical simulations individuated in common agreement that 
the EA Heisenberg model had only a zero temperature phase transi¬ 
tion, in evident contrast with experiments, which clearly showed a finite 
Tsg■ The scientific community used to explain this inconsistency in¬ 
voking the weak magnetic anisotropy present in real spin glasses, that 
wasn’t reproduced in numerical simulations, that caused a shift to the 
Ising behavior of the real materials. Nevertheless, this explanation wasn’t 
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at all satisfactory, since the critical exponents found in the experiments 
were completely different from the Ising ones. In 1992 an elegant theory 
was proposed by H.Kawamura, who suggested that a finite-temperature 
phase transition took place in the chiral sector. This scenario kept also 
after more modern simulation stated that the model has Tsc >0. It 
consists in two steps: 

1. Spin-Chirality decoupling 

(for a completely isotropic system, as the simulated ones) 

2. Spin-Chirality recoupling 

(for a weakly anisotropic system, as the experimental ones) 

The first part claims that the ordering is not simultaneous for both the 
order parameters (chiral and spin glass overlap). When we decrease tem¬ 
perature from the paramagnetic phase we first meet at Tqg the sponta¬ 
neous breaking of the discrete Z 2 symmetry with the preservation of the 
continuous SO (3) symmetry, and only at lower temperature T$g < Tea 
the whole 0(3) symmetry is broken. The higher transition, in which only 
the Z 2 symmetry is broken is called the “chiral glass transition”, and the 
phase between T C g and T S g, in which only the chiralities are ordered, is 
called the “chiral glass state”. 

In the second part of this scenario this result in adapted to a real 
spin glass, with inevitable weak magnetic anisotropy, which reduces the 
Hamiltonian’s symmetry from 50(3) x Z 2 to only Z 2 , mixing the chirality 
to the spin sector. So of the two transitions that take place in an isotropic 
system, it would be the chiral glass transition to dominate the real spin 
glasses’ phase transition. 

The chirality decoupling scenario explains in a natural way questions 
on the experimental phase transition like the origin of the non-Ising crit¬ 
ical exponents experimentally observed in canonical spin glasses, and it 
gave a very elegant explanation of what could be occurring when the early 
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numerical simulations estimated Tsg = 0. Yet, specially now that it is ac¬ 
cepted that Tsg = 0 was a problem of finite size system, or interpolations 
too far from the critical temperature, it is not a dominant scenario in the 
scientific community, since although there are recent works that would 
confirm it [16], there are others that give evidence of a simultaneous spin 
and chirality transition [17]. 

For what concerns the type of ordering, while for the SG order pa¬ 
rameter the nature of the glassy phase is debated not differently than 
with Ising spins, in the chiral sector there has been some evidence (figure 
2.2) of a one-step RSB [7], so the P(q) in the thermodynamic limit would 
be represented by two delta-function peaks at ±Qea and one at q = 0. 



Figure 2.2: The overlap distribution at T = 0.15 for the chirality of the 3D Heisen¬ 
berg spin glass with ±J couplings from [7, Kawamura (2005)]. This plot 
suggests that there may be a one-step RSB for the chirality. 


Mean Field Analysis Evaluating the help that the mean held 
resolution gives to the understanding of the finite-range interaction model 
is a troublesome issue. Table 2.1, which compares the mean held critical 
temperature T|^. F = —, where z is the coordination number (z = 6 in our 
case), with the estimates of the finite-range Tsg, shows us how increasing 
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m it becomes always less reliable. What happens is that the fluctuations 
neglected with the mean field approximation are more important the 
larger m is, so the mean field results with Heisenberg spins are even less 
reliable than those obtained with Ising spins. We remark also how the 
data in the table suggests that Tsg = 0 for m —» oo. Unfortunately no 
mean held solution exists in the chiral sector because it is impossible to 
define chirality in a non metric space. 


m 

Model 

rpMF 

1 SG 

Tsg 

T /rpMF 

2 SG/ ^ SG 

1 

Ising 

2.45 

0.97 [18] 

0.40 

2 

XY 

1.22 

0.34 [19] 

0.28 

3 

Heisenberg 

0.82 

0.12 [17] 

0.15 


Table 2.1: Comparison of the critical temperature, between mean field approximation 
and numerical evaluation, varying the number of components of the order 
parameter. We see that Tsa/T^f is small and decreases with to. Physi¬ 
cally this means that fluctuations get larger increasing to. Data suggests 
that Tsg —■> 0, and evidence has been found for this conjecture being 
true [20]. 

2.2 Dynamics 

At low temperature the dynamics in spin glasses become very slow, 
such that below the critical temperature the system is always out of 
equilibrium. There are many out of equilibrium effects in the dynamics of 
spin glasses that still haven’t been well understood as aging, rejuvenation 
and temperature cycles. Yet, in the present work we want to concentrate 
in another aspect that attracted our curiosity, that is the possibility of 
a straightforward analogy between the dynamics of spin and structural 
glasses. 

2.2.1 Energy Landscape Scenario 

Although out of equilibrium dynamics is proper of the spin glass 
phase, an interesting feature could be revealed if we notice that works 
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[8] that look for the critical temperature with dynamical criteria tend to 
overestimate it. This could suggest an analogy with structural glasses, 
in which Mode Coupling Theory predicts a power law divergence of the 
correlation time at T = T M q > T c . This divergence doesn’t take place 
because when we get close to Tmc different effects come in play, but it 
could be what is being observed in those works. Goldstein’s scenario 
[21], which by now is widely accepted, says that the dynamics at low 
temperature are completely different from those at high temperature. 

When T is low a supercooled liquid explores the phase space mainly by 
activated jumps between the local minima of the free energy. The char¬ 
acteristic time r to jump from one minimum to another is, in agreement 
with the Arrhenius law r oc , proportional to the exponential 

of the energy gap. On the contrary at higher temperatures the typical 
configurations are not local minima with a positive definite hessian, but 
saddles with negative directions, so the dynamics is much faster. The 
characteristic time does not represent any more how much it takes for 
the thermal agitation to cause the barrier jump, but how much it takes 
to the system to individuate the negative directions in the hessian. We 
understand, in this scenario, how the decrease of negative directions in 
lowering the temperature may increase the relaxation time making it 
apparently diverge, just until it becomes of the order of r oc exp -/3ziv . 

Our main aim in this work was to find evidence of the influence of 
the energy landscape in the dynamics of a spin glass, to give a basis 
for constructing a solid connection with the theory of the dynamics of 
structural glasses. In this perspective we did a deep investigation of the 
magnitudes that concern dynamics, such as time correlation functions, 
and tried to connect them to the energy landscape, by making a massive 
research of inherent structures (local energy minima). 
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3.1 Thermalization Dynamics 


3.1 Thermalization Dynamics 

Monte Carlo (MC) is the most used technique for spin glass com¬ 
puter simulations. It is a stochastic Markov Chain method that aims not 
to sample the whole configuration space, but only its most significant 
parts (importance sampling) in order to extract a fine estimate of the 
observables. 

For Ising spins there exist numerous MC variants that one can easily 
obtain imposing detailed balance to be satisfied, such as [22] the most 
famous Metropolis-Hastings, local sampling methods, perfect sampling 
and heat bath. A1 those consist in flipping sequentially the spins with 
the constraint of preserving the detailed balance. 

For multi-dimensional spins the deal is similar. We will now show the 
extension of a couple of those methods to Heisenberg spins. 


3.1.1 Metropolis 

The fingerprint of a Metropolis updating scheme is that new con¬ 
figurations are accepted with a probability min[l,e~ l3AE ]. We generate 
new configurations by modifying the orientation of a single spin in the 
following way: 

A3 — D random vector r in the sphere of radius 5 is extracted. 


The new spin proposal is 


^new 


sf + r 


~>old 


+ r1| 


(3.1) 


The new spin is accepted with probability P = min[ 1, e @ AE ] and re¬ 
jected with probability 1 — P. 


We see that the probability of proposing s~ d —> Sg ew is the same of 
proposing s^ ew —> s°^ d . The choice of 5 must be sufficiently accurate to 
give an acceptance rate of 40% ~ 60%, because a too small 5 (hence a 
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too small rejection rate) would make the dynamics too slow, while a too 
large rejection rate would keep the system still. 


3.1.2 Heat Bath 


The Heat Bath (HB) updating scheme has the main advantage of 
reducing to zero the rejection rate, and for this reason for Heisenberg 
spins it is much more performing than Metropolis. The main idea is 
that, chosen a site, the new spin proposal does not depend by the old 
spin, but it is chosen according to its probability distribution. 

For a given site the local held is known, so the conditional probability 
for the spin’s orientation is 


dP(s\h ) = ds- 


e fiJ(s-h) 

f dseP J PP 


(3.2) 


We can make a wise reference change to coordinates (V, y', z') with z' // 
h to get trivial scalar products, and then pass to polar coordinates, with 
.s = (cos O' cos 0', sin O' cos 0', sin (J)' = z'): 


dP(s\h) = d(j)d0- 


e /3JhsmO g i n q 


f d(f)dOed Jhsine sinO 
This way we get a factorized probability distribution function 


d(j) , , e 


dP(s\ti) = — x dz' 


fijhz' 


(3.3) 


(3.4) 


2ir f dz'eP Jhz ' 
where h = ||/?.||. Since the angle 0 < 4>' < 2tt is uniformly distributed 
it’s easy to get its pseudo-random value. For z! we use u = e^ Jhz ' that is 
uniformly distributed in the interval e~ /3Jh < u < e l3Jh . 

At this point we have to come back to our original frame with mini¬ 
mum computational effort. To do so we remark that our spin is 


s — z' — + y/l — z' 2 R± (3.5) 

where R± is a unitary random vector, orthogonal to the local held h, 
that we can obtain by extracting a random unitary vector, peeling its 
component parallel to h, and renormalizing it. 
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3.1.3 Over Relaxation 

The present is a deterministic, very computationally-fast, microcanon- 
ical simulation scheme. Although its dynamics are not physical, it helps 
when we want to sweep quickly the configuration space. Performing 1 HB 
sweep every L Over Relaxation (OR) sweeps we generate in the lattice 
an iso-energetic spin wave after each energy change. Of course it cannot 
reproduce real dynamics, but it is very useful to generate thermalized 
configurations of the system. 

The algorithm proposes the maximum spin change that leaves the en¬ 
ergy invariant so the new spin will always be accepted (see the Metropolis 
acceptance rate in section 3.1.1). The new spin is so binded the condition 

P ew -h = s° ld -h (3.6) 


and we obtain the maximal change by reversing the spin’s component 
that is perpendicular to the local field h: 


gnew 


= 2 h 


hx * '-’a: 

! ^T 


- S.7 


(3.7) 


3.2 Minimizing the Energy 


Obtaining spin glass low energy configurations is computationally 
hard [23], and this makes our objective of finding inherent structures a 
real challenge. While for Ising spins there exists a whole literature of 
numerical methods used to minimize the energy such as simulated an¬ 
nealing [24], cluster algorithms [25], genetic algorithms[26, 27], and com¬ 
binations of those methods[28, 29], almost nothing exists for Heisenberg 
spin glasses. 

We will now give a review of the small amount of methods used for 
Heisenberg spin glasses, and present our choice, which is not applicable 
to Ising spins and is new for HSG. 
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3.2.1 Gauss-Seidel 


The most trivial and common way to find low-energy states is per¬ 
forming a quench. 

It is easy to see that the energy of a single site 




—sg • h* 


(3.8) 


is minimized if the spin s x is parallel to the local field h x . this observation 
can suggest us a quite brutal way to lower the energy, i.e. aligning 
sequentially the spins to the local field, i.e. making the substitution 


s-r —y 



(3.9) 


This type of method is much older than the mere existence of spin glass 
models. In fact we will see that it is exactly the application to our model 
of the Gauss-Seidel method, used to solve systems of linear equations[30]. 
In fact the problem solving the linear system Ap = /, where A is a given 
symmetric positive-definite matrix and / is a given vector, is equivalent 
to the one of minimizing the Hamiltonian 


H(<f) = - (/>•#>)- 


(3.10) 


Given the Hamiltonian in equation 3.10, the Gauss-Seidel Algorithm con¬ 
sists in sweeping through the sites i, and at each step replacing pi by that 
new value (pi which minimizes H when the other fields {pj}j^i are held 
fixed at their current values[ 30]. In other words, if we write 

H((pi, (pjj fii) = ^Aifipi - bi)' 2 + a (3.11) 


the Gauss-Seidel algorithm consists in simply replacing pi with bi. This is 

exactly what we do with the quench routine when we substitute sg with 

ttA—. so if we show that the EA-Heisenberg Hamiltonian is reducible to 
IIMI 

the same form, we can state immediately this algorithm’s main property: 
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it’s ineffectiveness, and seek for a better solution. As a matter of fact this 
algorithm has an extremely slow rate of convergence [31]. For example 
for A = A + m 2 the convergence time grows as m~ 2 . 


Reduction of the EA Heisenberg Hamiltonian to a quadratic 
form Let’s show how, when we are close to a local minimum, our Hamil¬ 
tonian can be reduced to a quadratic form in chiral perturbation theory. 

If we call {s°} a generic inherent structure configuration, i.e. with all 
the spins aligned to the local field, any spin is expressible as 

sg = s°yjl -7?! + ^ (3.12) 


with the condition 

hx ■ TTx = s°x ■ Tz = 0. (3.13) 

7? is the component of the spin perpendicular to the local field, and we 
call it pion. If 17?| is small, i.e. we are close to the local minimum, we 
can develop the Hamiltonian to the second order in n, keeping in mind 
condition 3.13, in order to obtain 

H ({s}) = + JxgKy). (3.14) 

~ X x p-j7||=l 

This way the variation AH can be written as a quadratic form in an 
TV-dimensional space: 

AH = - +o(J 3 ) (3.15) 


with 


±D 

Aixy &xy ^ ( ^yx+e^Jxyi (3.16) 

/i==tl 


so finding an inherent structure is actually equivalent to minimizing a 
quadratic form when we’re in the limit of small pions. 


34 



3.2 Minimizing the Energy 


3.2.2 Genetic Algorithms 

Other than brutally applying the Gauss-Seidel method, very little has 
been done for Heisenberg spin glasses. Nevertheless, in a recent article 
[32] from Y.Iyama and F.Matsubara we can see an attempt to extend the 
genetic algorithm thought up by Pal [28] for Ising spins, to Heisenberg 
spins, with the intention of hireling ground states. Since this is the only 
other kind of energy-minimizing routine we found in literature we will 
explain the ideas on which it is based. 

The basic idea of a genetic algorithm is to find a low energy state 
mixing in a smart way the parts of a selection parent configurations 
(from analogy with reproduction the name genetic, since it is as if the 
configurations pass to their heirs their genetic code). Clearly the initial 
choice of the parent configurations is a crucial step. 

Pal’s suggestion [28] is to use a hybrid algorithm which alternates 
the genetic algorithm steps with a generic non-genetic minimization al¬ 
gorithm, and this is what is done by Iyama and Matsubara in [32], who 
chose to alternate quenches (Gauss-Seidel) to the genetic steps. 

It appears, from their work, that a very big computational effort is 
necessary to find low energy configurations with their algorithm, since 
they study very small lattices (L < 13), and the number of samples 
studied is not high (500 samples for L = 13). 

Later in this paper we will refer to this work since it is the only ex¬ 
ample in literature of search of ground states in Heisenberg spin glasses 
(although they study ±J bonds, while ours have a Gaussian distribu¬ 
tion), and they give the only existing (to now) estimate of the ground 
state energies. In a brief study of ±1 bonds we will show how results 
comparable to theirs are achievable with much smaller computational 
effort. 
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3.2.3 Successive Over Relaxation 


It is common sense that greedy algorithms such as Gauss-Seidel are 
not very effective, especially compared to ungreedier versions of their 
selves. Taking advantage of this precious advice we can improve this 
scheme in a simple and effective way. Instead of choosing the fa that 
minimizes the Hamiltonian we can make an apparent non-optimal choice 
that lowers the energy, but less. On the long run the convergence time 
is definitely quicker (for A — A + m 2 the convergence time grows as m~ 1 
instead of m~ 2 with the best choice of the parameters). 

In our problem, if we call 

r -gs _ hg 
3 II hg 

the new spin that minimizes the local Hamiltonian, and 


(3.17) 


rOH 

a x 


v hg 

= 2 h. 


hi 


(3.18) 


the spin deriving from the Over Relaxation algorithm described in section 
3.1.3, the Successive Over Relaxation is expressed by choosing the new 
spin as 


jSOR _ 
* \ Z — 


s9 s + A s9 R 


’ A '* ||jf s + Asf s || (3 ' 19) 

where A > 0 is the parameter that determines how ungreedy the rou¬ 
tine is. This way we are choosing our new spin damping the maximum 
energy minimization with the iso-energetic change that most varies the 
Hamiltonian. 

Figure 3.1 shows how the energy decrease is initially lower with higher 
A’s, but there is a time after which the ungreedier routine (greater A) 
becomes more efficient. 

We remark the two trivial limits in A for this algorithm: 


lim s' 

A—aO 


SOR 

A 




lim 


A —yoo 


= s'- 


(3.20) 

(3.21) 
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Evolution of the energy in a L=100 lattice with t =0 



Figure 3.1: Successive Over Relaxation in a L = 100 after having performed t* = 0 
Gauss-Seidel steps before starting with SOR. With small A the dynamics 
is similar to a quench: the system dives into a valley of the free energy 
without selection, and it stays there a lot of time before converging, as 
we know that with the Gauss-Seidel method convergence is a bad issue. 
On the other side dynamics with very large A spend a lot of time at high 
energy: they select better the valley where to get inside, so they reach 
lower energies (it will be more evident in the zoomed figure 4.2) and once 
they’re in a valley they converge faster. We can remark the too greedy 
routines, which converge too slowly because they spend too much time at 
low energies, and the too greedy ones, that stay too much time at high 
energies. 

It is clear that A has to be chosen with care, because a small one risks to 
be too greedy, hence too similar to a brutal quench, while a too big one 
may be too ungreedy, having a too long plateau, and never converge. 
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Variants If we want our routine to be more site-dependent we can 
slightly modify it as follows: 


^new 

^X,x 


hg + A Sg R 
\\hg + Xs9 R \\ 


(3.22) 


This way the dynamics will be less damped in the sites where the local 
field is strong, and more where it is weak. This is the variant we used 
throughout our investigation. 

Another idea could be to alternate aL (a = 1,2) Over Relaxation 
steps with one quench steps instead of updating the spins with always in 
the same ungreedy way. This would mean that before lowering the energy 
we sweep the system with an iso-energetic wave of spins. We tested this 
method both alternating GS and OR, and SOR and OR. The result 
was more satisfying than pure Gauss-Seidel, but the convergence time 
resulted being less convenient with respect to Successive Over Relaxation. 


3.3 Computing Observables 

A main advantage of simulations is that we can know the exactly 
properties of each single particle of the system, while in real experiments 
only some (global) observables can be found by analysing the reaction to 
determinate stimuli. The possibility of computing observables directly 
represents a great advantage and it pays off the problem of small sizes. 

The observables that we computed in our simulations were: 

• Energy 

• Overlaps 

• Correlation Lengths and Functions 

• Magnetic Susceptibility 

• Energy-Overlap Correlation 
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• Binder Cumulants 

• Correlation Times and Functions 
Let’s analyse them one by one. 


3.3.1 Energy 


The energy is the most intuitive and basic observable to look for in a 
simulation, and it represents one of the great challenges of SG systems, 
since clue to the frustration it is very hard to hnd low energy configura¬ 
tions, and in fact this is what we’re going for. 

It can be found very trivially, as it is defined by the model’s Hamil¬ 
tonian: 

H = ^ Sg ■ h s (3.23) 

X 

It simply requires a sweep over all the spins summing the energy per site: 





(3.24) 


3.3.2 Overlaps 


There are two types of overlap to compute, the CG overlap, and the 
SG. We used the overlap definitions anticipated in chapter 1. The CG 
overlap was the scalar quantity 

3 3 

£<•(*), 0 . 25 ) 


while for the SG overlap we used a slightly more complicated object, 
also introduced in chapter 1, that guarantees rotational invariance. As a 
matter of fact it is not an overlap, but a square overlap: 


i- = AmW)- 


(3.26) 


39 



3.3 Computing Observables 


Properties of q 2 Let’s see the typical values of q 2 . We remind that 
the tensors Q a p are defined as 

Q«p = r aP {x) = s { £ a sf^ (3.27) 

X X 

so, by substituting this definition and rearranging the sums we get, with 
trivial calculation, 

tr(QQ') = E E(4“’ - 4 a) )(4‘ ) ■ 4’b < 3 ^ 28 ) 

x y 

that may be easily bounded thanks to the Cauchy-Schwarz inequality, 
that guarantees us that sgsp < ^Js = 1, so 

tr{QQ ] ) < IV 2 , (3.29) 

and the square overlap, with the normalization that we chose, gets the 
bound 

q 2 =< 1 (3.30) 

SG Self-overlap We want to show the typical value of the SG self- 
overlap for a random configuration. To evaluate it we impose (a) = (b ) 
in eq. 3.28, and obtain 

t' r (QselfQ\ e lf ) = E Eft ' %) 2 (3.3i) 

x y 

= E4+Eft'%) 2 - < 3 - 32 ) 

x Xy 

To estimate the average of the product (s# • s^j) 2 we can put ourselves in 
a reference frame (x, y , z) with sg/ / z. The value of the scalar product is 
then given by the third component of Sjj, s$ jZ . This way 

(sg ■ s$) 2 = sl z (3.33) 

is a uniformly distributed random variable whose mean value is 1/3, since 
our system is completely isotropic. 
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We can now determine the average self-overlap 

^ Qself ' > ~ 772 T ^ ^ < (Sf • Stf) >) 

x^y 

J_ , 1_ l_ 1 

N ^ 3 3 N * 3 


(3.34) 

(3.35) 


3.3.3 Correlation Lengths and Functions 


An extremely important quantity for our simulations was the correla¬ 
tion length, because not only it gives information regarding the properties 
of a state (e.g. it is zero at infinite temperature while it diverges at the 
phase transition), but since it estimates until which distance two spins 
can be considered correlated, it indicates us if we can consider the system 
in thermodynamic limit or not. For the latter check we can also see if 
the correlation functions go to zero within half lattice size 1 . 

The correlation function relative to a generic observable 0 is defined 
as: 

G(r) = ^ (px<px+? (3.36) 

X 

As it is shown in appendix A.l, together with a computationally smart 
way to write it, a way to compute the correlation length £ involves the 
correlation function’s Fourier transform G(k): 


Zf 2 = 


G( 0) 


Asin\lYG{k. n 


- 1 ] = 


4 sin 2 (f) L F 


1-1] 


(3.37) 


where we stressed that the magnetic susceptibility y is the k = 0 Fourier 
transform of the correlation function: 


X = G(0) 
F = G{k min ) 


(3.38) 

(3.39) 


and due to the periodic boundary conditions k min = (^4,0,0). 

1 Since the boundary conditions are periodic the maximum distance between two 
points is Lj 2, not L. 
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Integral Estimators If we define an integral estimator //. as: 

L 

I k = f 2 dfr k G(r) (3.40) 

Jo 

we can readily see that (appendix A.2) in the thermodynamic limit we can 
define two more ways to find the correlation length in a D-dimensional 
system: 

£<° 2 ) ( 341) 

y 1 d -i 

ZkM 1 = ^ cx e (3.42) 

All these definitions suffer from systematic errors because they are 
derived from the long distance decay of G: 

G(r)~4 e *P(-;) < 3 - 43 ) 

that is only an asymptotic formula for large values of r. The systematic 
errors in these definitions can be reduced by considering a large value of 
k (since the r k factor would suppress the deviations at short distances). 
However, there is also the issue of statistical errors to consider. As it 
can be seen in figure 3.2, borrowed from reference [33], a large value of 
k pushes the maximum of r k G{r) towards a zone where the statistical 
error becomes too big. Therefore it is necessary to find a compromise k. 


Plane to plane correlation functions We can define three types of 
plane to plane correlation functions, one for each direction. To avoid 
making notation heavy we write it explicitly for a given direction: 

1 L_1 

C x (r) = — P(x)P{x + r) (3.44) 

x=0 


-P(x) ^ ^ 4*xyz 


(3.45) 
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Figure 3.2: Comparison between r 2 G(r) and r 4 G(r) from simulations of an L = 80 
Ising spin lattices at three different temperatures. In the first case the 
significant data is clearly out of the noise-dominated region, while in the 
second the maximum is pushed rightwards giving plenty of noise to the 
calculation of J 4 . Figure borrowed from [33]. 


where (j> xyz is any observable with which a correlation function has sense. 
Since the space is perfectly isotropic we can reduce the error by calculat¬ 
ing 

C(r) = \(C x (r) + C y (r ) + C z (r)). (3.46) 

o 

Plane to plane correlation functions represent a fine way to have more 
precise estimates of the correlation length £, since beyond being a nat¬ 
ural choice for square lattices, in which we can easily average it on the 
euclidean directions, and they decay far more slowly than the point to 
point correlation functions. 

I 11 fact if we define the plane to plane integral estimators: 

J k = 1 2 drr k C(r ) (3.47) 

Jo 

where we called C(r) the mean of the plane to plane correlation func¬ 
tions in all the directions, it is immediate to see that since the space is 
homogeneous and isotropic 


Jk ~ 47r/ fc+2 
X = 2 J 0 - C(0) - C(L/ 2) 
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(3.48) 

(3.49) 
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so for the same integral estimator we need to use lower powers of r if 
we’re using the Jj-’s. This means that r k C(r ) will have sense in a larger 
span of distances, so the amount of data without signal will drastically 
decrease. 

Truncating data Another thing that can be done to reduce con¬ 
sistently the statistical error, in favor of a small bias, is to do as proposed 
by Mart in-Mayor, Marinari et al. in [33], i.e. to use a self-consistent in¬ 
tegration cut-off as it has always been done with the study of correlated 
time series [34]. 

What we do is to decide a posteriori a critical distance r c beyond 
which collecting data has no sense because the measure is dominated by 
the statistical error. This distance was chosen to be the first one at which 
the error became greater than one third of the correlation function: 

L 

2 r c 

h = J2 ^ G(r ) —> 4 = E ^(0 ( 3 - 5 °) 

0 0 

r c = {first r : C(r + 1) < 3<r(r + 1)} (3.51) 

A truncation of this type in our work has reduced the error of several 
orders of magnitude, giving a quite small bias [33] (one or two percents), 
that could be seen from the change of the value of the magnetic suscep¬ 
tibility y. 

3.3.4 Magnetic Susceptibility 

In theoretical physics, the magnetic susceptibility is defined as the 
zero-moment Fourier transform of the correlation function 

x = C(k = 0) = Y / C(r) (3.52) 

r 

The magnetic susceptibility that we calculate in our simulations is 
what in the experimental jargon is called the non-linear susceptibility. If 
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we put a spin glass under a uniform magnetic field h, the magnetization 
m can be written as a series of h 


m{h) = XLh + XNLh 3 + ... (3.53) 

The term \l is the linear magnetic susceptibility ^ that does not diverge 
at T S g ■ What diverges is the non-linear magnetic susceptibility \nl and 
that’s what makes it a more interesting observable. 


3.3.5 Binder’s Cumulants 

Binder’s cumulants are ratios of cumulants of the order parameter 
that give us a good comprehension of the finite-size effects. They are 
often used to individuate phase transitions because at T c they are size- 
invariant. In our case we used them to have a check of how close to the 
thermodynamic limit we were throughout our work. 

Due to the presence of two order parameter, we can define two Binder 
Cumulants. The chiral-glass Binder Cumulant is 


B cc = < E ca> 


< E ia> 


(3.54) 


while the spin glass cumulant is 


B? g = 


< tr(QQ*) > s 


< q 2 > 2 

< q 4 > 


(3.55) 


< br(QQt) 2 > 

It is possible to show that the asymptotic values of these cumulants 
for systems at infinite temperature in the thermodynamical limit are 

. B GG = 11/9 

• B£ G = 3 

The interested reader can refer to appendix A.3 for the calculation of 
B sg . 

OO 
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3.3.6 Self-correlation Times and Functions 

The self-correlation time r represents the amount of time in which a 
system loses all the memory of its initial configuration. It is calculated 
from the time-correlation functions C(t), that give us an idea of how 
much the system at time t is correlated with itself at time 0. The time- 
correlation functions measure how much a particular observable changes 
during time, so it is evident that it s possible to define a whole set of 
self-correlation times and functions, and that the most important of the 
times is the greatest, since it is the characteristic time after which all the 
observables have no correlation with their initial values. 

In our work these quantities have been used both to know how much 
we needed to wait for the thermalization of the system, and to measure 
the scaling of the r in the paramagnetic phase. 

We used three self-correlation functions: 

• Scalar Correlation function 

CScalar (t) =< Sj?(0) • Sg(t) >= ^ Sj?(0) • Sg(t) (3.56) 

X 

• Tensorial Correlation function 

C'Tensorialit ) 2 ^ (Qo,tQoj ) (3.57) 

where Qo,t is the overlap calculated between configurations of the 
same replica at times 0 and t. 

• CG Correlation function 

C*Scalar(t) y £ £ 4(0)<-« (3-58) 

fl X 

And calculated the correlation times: 
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T Scalar / CScalar(t)dt (3.59) 

Jo 

/»oo 

t C g = / Ccdfydt (3.60) 

Jo 

Just as in the case of the spatial correlation times, we truncated the data 
when C(t) was smaller that 3 times its error bar. We did not calculate 
TTensoriai because the square overlap does not go to zero for uncorrelated 
samples, but to 1/N, so the data would be biased and not L-independent, 
especially when truncating data. 

Aging It looks obvious that if one measures the self-correlation func¬ 
tion starting not from an initial time 0, but from a time to, nothing should 
change. Nevertheless, in spin glasses this time-translational invariance is 
not so obvious. In fact what happens is that while this symmetry is 
verified at high temperatures, below Tgc the value of the self-correlation 
function depends on the amount of time we wait before starting our 
measures. This phenomenon, called aging, obliges us to redefine the cor¬ 
relation functions as C(t w , t ), where t w is the time we wait before starting 
to collecting data. 

We measured self-correlation functions from many t w , to verify that 
we did not have aging. Indeed the data was independent from t w as 
we can see from the collapse of the data in figures 3.3 where we plotted 
C(t] t w ) for many t w . It is because C(t w , t ) = C(t — t w ) that in this report 
we refer to it as C(t). 
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Absence of aging of C sca | ar at T=0.3 for L=64 



Absence of aging for C CG at T=0.3 for L=64 
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Figure 3.3: The scalar and chiral-glass time-correlation function C sca i ar (t w ,t) and 
CcG{tw,t) in an L = 64 lattice at T = 0.3. As we can see there is no 
aging in the paramagnetic phase, since data for different t w ’s collapse 
perfectly. 


48 







Chapter 

Results 


Contents 

4.1 Using Gauss-Seidel . 50 

4.2 Using Successive Over Relaxation. 51 

4.2.1 Parameter choice. 51 

4.3 A comparison with genetic algorithms .... 56 

4.3.1 Inherent structures at infinite temperature . . 58 

4.4 Inherent structures at finite temperature . . 62 

4.5 Dynamics. 71 

4.6 Extrapolations. 72 

4.6.1 An Estimate of Tsg . 73 

4.6.2 An Estimate of the Exponent v . 74 

4.6.3 Self-correlation times. 75 

4.6.4 A Lower Bound for the Dynamic Exponent z 77 

4.7 Minima or Saddles?. 78 


49 













4.1 Using Gauss-Seidel 


4.1 Using Gauss-Seidel 

Although we showed in chapter 3.2.1 that the Gauss-Seidel algorithm 
is not effective for our problem, a complete analysis should start from the 
simplest cases to be able to make comparisons which make us understand 
if different routines really do improve the dynamics or not. 

In this perspective we investigated the convergence of the Gauss- 
Seidel algorithm in quite large lattices, and as it is shown in figure 4.1 
the convergence rate is with a power law, which is far too slow for our 
needs. 


HSG 3D L=100 quenches from random configs 



Figure 4.1: Evolution of the energy density during a quench in a L = 100 lattice. The 
plot shows 20 quenches from random configurations (red lines), i.e. from 
T = oo configurations, which were stopped when the energy gain after a 
sweep was less than 10 -10 . The green line is the mean with error bars, 
which was performed only in the points where all the data was present, 
since the duration of the runs was not even. The dashed blue line is a fit 
in the linear regime f -1 < 0.0005, which yields a good y 2 test. 


At the beginning of our work our convergence criterion has been 
checking the energy gain in the last sweep, but we readily realized that 
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a more solid one was to see the maximum pion in the lattice, since in an 
inherent structure the pion held is zero. 

4.2 Using Successive Over Relaxation 

A disadvantage of Successive Over Relaxation respect to classical 
quenches is that it depends from the parameter A, therefore it’s per¬ 
formance is not univocally determined, and it is necessary to do some 
preliminary studies to choose the optimal one. We demand a good algo¬ 
rithm the capability of converging to low energies in an affordable amount 
of time. 

Since at the beginning of a simulation a greedy routine such as Gauss- 
Seidel descends in energy much faster than an ungreedy one such as Suc¬ 
cessive Over Relaxation, one could argue that to increase the convergence 
speed it could be convenient to start with the Successive Over Relaxation 
only after having performed t* quench steps. This increases the number 
of non-physical parameters of our routine to two: A and t*. 

4.2.1 Parameter choice 

To decide how to choose t* and A we ran simulations on a large lattice 
seeking for the best performances. 

We start by stressing the properties of t*. Figure 4.2 shows how the 
convergence time does not depend from this variable, but completely on 
A. We can also see (figure 4.3) that, as intuition suggests, more Gauss- 
Seidcl steps we perform before starting with A / 0 less we will be able 
to descend in energy: it is impossible for the Successive Over Relaxation 
to recover the “unwise” greedy quench steps. 

So if making a number of quench steps before does not decrease the 
convergence time, and it entails an increase in the final energy, there is 
absolutely no point in keeping t* ^ 0. 
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Dependence of e* IS) on X and t in a L=100 lattice 



Figure 4.2: The evolution of the energy for 3 different A’s, and t*’s in a L = 100 lattice. 

Each A has an associated color. The runs stopped when the maximum 
spin component perpendicular to the local field was smaller than 10 -10 . 
It is possible to see the point t* = 100 where the Gauss-Seidel routine 
stops in favor of the SOR. Depending on A, after this point, the evolution 
at this point is faster or slower than the quench. The most ungreedy 
A on the long run reaches a lower energy, but at the cost of a much 
longer convergence time. What appears clear is that the convergence 
time depends exclusively on A, and also that the greater the quantity of 
quench steps before turning to SOR, the higher are the energies achieved. 
Hence, there is no point on having t* > 0. 


The situation is more delicate when we proceed to select an optimal 
A. In this case the choice is not univocal, since it depends on which 
criteria we consider dominant. 

Looking at the raw data one could be attracted by the possibility of 
a pretty scaling law for A such as, for example, 

e.( IS \\) = + A\~ a ', (4.1) 

from where it would be easy to extrapolate the energy of the inherent 
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Energy density vs t* _1 



Figure 4.3: The final energies one obtains with different t*’s in a L = 100 lattice. 

Each color indicates a different A. Definitely the energy of the inherent 
structure is higher the greater t* is, no matter A. 

structure. Unfortunately it is not so direct to state the existence of a 
similar scaling law (fig 4.4), and also if it exists, it is not as simple as 
the one described in eq. 4.1. Since the objective of this work is not to 
find all the characteristics of our minimizing routine we did not linger in 
finding a scaling law. 

The other criterion to keep in mind when choosing our A is the con¬ 
vergence time, which is absolutely not negligible. In fact although for 
small A the Successive Over Relaxation makes the convergence faster, 
there is a A around 100 (see fig. 4.5) in which the algorithm becomes too 
ungreedy, and the convergence times start growing. Indeed it is true that 
any A < 10 5 represents an objective gain respect to the quench, since it 
yields lower energies in less time. 
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Figure 4.4: The energy of the inherent structures in a L = 100 lattice (except the 
point at A = 300 that comes from an L = 64 lattice). We can ideally 
divide the plot in two sectors. In the left sector we have large A’s: the 
dynamics takes place mainly at high energies and the system takes its 
time to choose the valley in which it wants to lay down. In the right 
sector, for smaller A, the system dives quickly in a minimum, losing the 
choice of a lower valley. Mind that the dependency of e l ' IS ' > on the system 
size L is studied in Fig. 4.7, below. With our statistical accuracy, the 
results seem L -independent from L = 32. 


Nevertheless we ought to choose only one single A between the infin¬ 
ity available. Since we did not believe that our algorithm could cross 
the phase transition for a specifically big A 1 , and start finding sub- 
exponentially numerous low energy configurations, and since the defini¬ 
tion of inherent structure is intrinsically related to the minimizing routine 


1 If the algorithm were to cross the phase transition, the correlation lengths would 
depend strictly on the system size for any L, since at the phase transition £ becomes 
infinite. On the contrary, we see that the correlation functions collapse perfectly once 
the system is large enough (figures 4.8 and 4.13), giving an L-independent correlation 
length £. 
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100000 
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Figure 4.5: Convergence times for different A. Here, as in figure 4.4, it is possible 
to identify the two regimes. For small A the system drops quickly into a 
minimum, but it takes it a long time to converge once in the valley. For 
large A it stays plenty of time at high energies. We see that from the pure 
Gauss-Seidel the convergence time decreases linearly with A, because it 
becomes more and more effective once in a valley. On the other side, if 
we increase A too much the algorithm becomes too ungreedy, loses too 
much time at high energies, and the convergence time starts growing with 
A. All the points come from L = 100 systems, except A = 300 that is 
extracted from L = 64. 


used, we had no problem choosing the A which made our simulations the 
fastest, for the rest of the investigation, except in the next section, where 
we wanted to remark the differences (other than the energy) between 
configurations achieved with different A’s or where we explicitly needed 
to achieve a very low energy. Another good reason for having chosen 
A = 100 is that a greater one would have yielded inherent structures 
with too long correlation lengths, preventing us from being able to calcu- 
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late inherent structures at the finite temperature we were interested in. 
For those reasons we made our investigation with A = 100, comparing 
it with A = 300 to see if there is some physical difference between the 
regime in which the system passes most of the time at low energies and 
the one where it delays at high energies. 

4.3 A comparison with genetic algorithms 

A nice check for our new algorithm is comparison with different meth¬ 
ods for finding low energy configurations. The only work of this type that 
we found in literature is [32], which uses a hybrid genetic algorithm to 
find ground states in EA Heisenberg spin glasses with ±J couplings. In 
spite of those two differences, we thought that it might be interesting 
to make a comparison with their results, so we made a little parenthesis 
in our work and studied the ±J model, since the maximum lattice size 
analysed by them has been L = 13 and this allowed us to not spend too 
much CPU time in this intent. 

We made a series of simulations with growing A to see how close we 
could get to what they claimed to be the ground state energy (we will call 
it e^ 4 ' 1 ). As the reader can see in figure 4.6 although we never reached 
i the gap between our lowest energy and theirs is of the order of a 
standard deviation. 

We also tried to reach lower temperatures by thermalizing the sys¬ 
tem at low temperatures and starting SOR from these configurations. 

Although they were slightly lower, they still were not in the error bar of 

JGA). 

e GS ■ 

1. e ±J (T = 0.14, A = 10 5 ; L = 13) = -2.0394 ± 0.0002 

2. e±J (T = 0.12, A = 10 5 ; L = 13) = -2.0396 ± 0.0001 

3. e ^ s A) = -2.0403 ± 0.0003 
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Figure 4.6: Inherent structures in the ±J Heisenberg model for L = 13. The energy 
density is normalized in conformity with ref. 4.6, i.e. without dividing 
by the number of dimensions. The blue points represent the average 
energies of the inherent structures we found for each A. The continuous 
black line represents obtained in ref. 4.6 with a hybrid genetic 

algorithm, and the dotted lines its uncertainty. The difference between 
the energies achieved by their protocol and ours is fairly small, and they 
can’t cope with system sizes greater than L = 13, where finite-size effects 
are dramatic, so to our scope their algorithm would not be as useful as 
SOR, that converges quickly and permits us to study very large lattice, 
keeping our measures in the thermodynamic limit. 


It is to be noted that the correlation lengths of the inherent structures 
found after having thermalized at T < oo were appreciably different from 
those from T — oo. This means that the situation is not so simple as it 
may appear, and a state cannot be characterized by the only energy. 

We know that our relaxation routine is definitely not going under 
the phase transition, since the correlation lengths stay constant when we 
increase the size of the lattice. 
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If is really the ground state energy this means that the energy 
gap between the topological transition and the ground states is very 
small and the topological transition does not have the importance one 
would like to give it. Nevertheless this would demonstrate that the energy 
landscape is quite trivial until very close to the ground state, and a simple 
relaxation algorithm as SOR is enough to almost reach it. 

It could also be, on the contrary, that does not indicate the 

ground state energy, but it is only a different estimate of the energy at 
the topological transition. In this scenario the energy landscape would 
be quite non-trivial, proving that neither a hybrid genetic algorithm can 
easily pass the topological transition. 

Unfortunately in [32] the only physical magnitude reported is the en¬ 
ergy, so this does not allow us to get to any conclusion. We can only 
say that SOR has revealed quite satisfying, since it yielded, with a com¬ 
petitive CPU effort, energies at two standard deviations from what are 
claimed to be ground states. 

4.3.1 Inherent structures at infinite temperature 

The first physically relevant survey to do is the investigation of the in¬ 
herent structures of the typical states of infinite temperature, doing SOR 
from random configurations, since at T — oo all the configurations are 
equally probable. We stress immediately that, no matter the system size, 
for different A’s not only the energy of the states is distinct (figure 4.7), 
but also the correlation functions, consequently the correlation lengths 
and the magnetic susceptibilities, are different. This follows physical in¬ 
tuition, since a lower achieved energy corresponds to greater ordering. It 
is then clear how the concept of inherent structure is strictly binded to 
the routine we use to descend in energy. 

Table 4.1 shows the typical values of the correlation lengths and mag¬ 
netic susceptibilities one could extract from correlation functions such as 
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Figure 4.7: The energy of the inherent structures for different system sizes. Already 
from L = 16 this magnitude seems to become quite stable. Neverthe¬ 
less we’re far from the thermodynamic limit, since we’ll see that other 
quantities need bigger lattices. 


the ones in figure 4.8, for four different A’s. As we just argued, it is not 
surprising that the correlation lengths increase monotonically with A. 

Regarding the size of the system, from figure 4.7 we see that already 
from L — 16 the energy seems to be stable, so for T = oo we could 
consider ourselves in the thermodynamic limit from L = 32. Nevertheless 
we need to be careful, since energy is a quickly converging quantity. In 
fact from figure 4.8 we see that for small L’s the correlation functions do 
not go to zero. 

This fact is confirmed by the binder cumulants, which only for larger 
sizes are compatible with their asymptotic values, as it is shown in table 
4.2, where we compare inherent structures at infinite temperature for 
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Spin Glass Correlation function 



Figure 4.8: Normalized SG plane to plane correlation functions of the Inherent Struc¬ 
tures for many A’s and L 's at T = oo. The data was truncated when 
it became smaller than three times its error, since at that point it was 
dominated by the statistic noise. It appears clear from here that the 
configurations obtained depend strictly A since we can appreciate three 
different trends, one for each A. Moreover, this plot shows the strong 
finite-size effects for smaller lattices, since many correlation functions do 
not converge to zero. Notice how greater A’s require larger lattices. 


A = 100 for different system sizes. 
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Table 4.1: £ and % from T = oo, for different A’s. 


A 

isc 

icG 

Xsg 

XCG 

10 

100 

1000 

10000 

2.79 ±0.05 

3.06 ±0.04 

4.49 ± 0.07 

5.90 ±0.13 

0.93 ±0.03 

1.06 ±0.02 

1.57 ±0.05 
2.48 ±0.10 

178 ±3 
205.3 ± 1.4 
443 ±6 
927 ±22 

0.0992 ±0.0011 

0.1106 ±0.0010 

0.213 ±0.003 

0.432 ±0.021 




Table 4.1: Growth of the correlation lengths and the magnetic susceptibilities of the 
Inherent Structures from T = oo for increasing A. Higher A’s yield lower energies 
which correspond to more correlated configurations, as the data shows. All of the 
data comes from L = 128 lattices, except the one for A = 10000 that has L = 64. Of 
course it suffers fair finite-size effects, since we see from figure 4.8 that L = 64 isn’t 
enough also for A = 1000. 


Table 4.2: Binder Cumulants for A = 100 and T — oo. 


L 

b c l g 

B S L G 

8 

3.58 ±0.07 

1.071 ±0.002 

16 

3.34 ±0.09 

1.179 ±0.004 

32 

2.97 ±0.06 

1.215 ±0.004 

64 

3.05 ±0.06 

1.232 ±0.004 

128 

2.80 ±0.22 

1.213 ±0.015 

oo 

3 

1.2 




Table 4.2: The binder cumulants reflect the conclusions on the finite-size effects that 
we took when looking at the scaling of the energy with L. They tell us that from 
L = 32 we are already in the thermodynamic limit, since they are compatible with 
their infinite limit. The uncertainty for L = 128 is greater because we took less 
measures than with the other sizes, since there was no point in spending too much 
CPU time in data we did not really need for our investigation. 
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4.4 Inherent structures at finite tempera¬ 
ture 

To find inherent structures at finite temperature it was necessary to 
thermalize the system and only then go down in energy. Our ambition 
was to be always in the thermodynamic limit, so that thermalization 
could become a very time-expensive task for low temperatures. We chose 
temperatures from T = 0.19 and greater because a previous article based 
on dynamical studies [8, M.Picco and F.Ritort (2004)] affirmed that at 
that temperature a phase transition was taking place. Also, we wanted 
our investigation to be in the thermodynamic limit, so we could not 
descend too much in temperature because the closer to T S g we are, the 
longer the correlation lengths are. 

To get to know the thermalization time of the system we have waited 
until the energy became stationary and the three time-correlation func¬ 
tions described in section 3.3.6 went to zero. 

It turned out that although the thermalization time was not fast, 
it was not either prohibitive, so we have been able to thermalize 1000 
samples for each chosen temperature. For T > 0.23 a L = 64 lattice was 
enough, while for lower temperatures the correlations grew too much and 
it was imperative to use L = 128. 

Trend of the Inherent Structure Energy Figure 4.9 shows how 
e ( IS ) varies with the inverse temperature /3 = 1/T. It is almost constant 
from T = oo to T = 0.5, as the error bars of those two points overlap, and 
it starts descending appreciably around T = 0.3 (/3 = 3.3). This is the 
same qualitative scenario we encounter in structural glasses [35], where 
e <yIS ' > is temperature independent for high T’s, and it starts decreasing 
significantly only from some certain temperature. 
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Figure 4.9: Energies of the inherent structures in function of the inverse temperature 
0 = 1/T for two different A’s. We plotted in function of 0 to be able 
to show the infinite-temperature limit. The trend is qualitatively the 
same we encounter in supercooled liquids and structural glasses, where 
the energy of the inherent structure is almost independent from T at high 
temperatures, and it descends steeply at lower temperatures. Keep in 
mind that being in function of 0 the energy descent at low temperature 
looks less steep than what it would with T on the rr-axis. Nevertheless 
this lets us appreciate the extremely small variation between T = oo 
(0 = 0) and T = 0.50 {0 = 2), since their error bars cross, in contrast 
with the appreciable descent for lower temperatures (higher 0). Notice 
how the difference between the two protocols decreases with decreasing 
temperature, so the choice of A becomes less relevant the lower T is. The 
three points at lowest temperature (highest 0) have less uncertainty than 
the others because the low-temperature data requested using L = 128 
while for the other we used L = 64. 


Anharmonic Energy Term From the equipartition theorem we can 
divide the average energy at a given temperature in many contributions. 
There is a term given from the topology of the energy landscape, a har¬ 
monic term, proportional to T, that represents the thermal agitation 
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Figure 4.10: The anharmonic contribution to the energy. The red points refers to 
the left y-axis, and give the absolute value of the anharmonic term of 
the energy. The other lines refer to the right y-axis. The blue points 
indicate the relative weight of the anharmonic contribution. The value 
1 indicates that it has no weight, and that is why there is a horizontal 
line indicating it. From the blue points we understand that the anhar¬ 
monic term is no larger than the 12% at T = 0.5, which is way over the 
critical temperature. At lower temperatures, though always in the deep 
paramagnetic phase, the anharmonic contribution is less than half. This 
gives space to an interpretation of the energy of a state as principally an 
inherent structure plus harmonic thermal fluctuations. 


around a minimum, and an anharmonic term that permits the move¬ 
ment through different valleys, which is more relevant the higher T is. 
The relation is 

e(T) ~ e (/s) (T) + - + e {anh) (T ) (4.2) 

3 

where e ( ' anh \T) is the anharmonic contribution to the energy and the fac¬ 
tors 1/3 derive from the energy normalization we chose. We can compute 
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the anharmonic contribution to the total energy by reversing equation 
4.2, and compute 

e {anh )(T) = e(T) - e [IS) {T ) - - (4.3) 

3 


that alone does not mean a lot to us, so we can understand its relative 
weight by looking at the quantity 


e{T)-e {IS \T) _ ft («"'0(y) 

773 “ + T/3 


(4.4) 


Figure 4.10 shows the absolute anharmonic term of the energy on 
the yl-axis, and it shows its relative weight on the 2 / 2 -axis, where 2/2 = 1 
marks the absence of an anharmonic contribution. We see that its weight 
does not exceed 12% even when the temperature is quite high, since 
T = 0.5 ~ 5 Tsg, hence it represents a non-primary effect. This fact 
lets us figure out that the energy of a configuration is almost completely 
determined by the shape of the energy landscape plus harmonic ther¬ 
mal fluctuations over it. The small , but yet absolutely not negligible, 
anharmonic effects are those that make us find different energies when 
we change the algorithm of energy descent, since they’re the cause of its 
movement from one valley to the other. 


Correlation lengths We also want to show (figures 4.11 and 4.12) how 
correlation lengths vary with temperature. From these trends we made 
extrapolations that we will show in a further section. By now we can 
appreciate the apparent divergence that both present around T ~ 0.15. 

A more interesting fact that the correlation lengths can remark to us 
is that more we go down in temperature, more important the role of the 
energy landscape becomes. In fact we can see in figure 4.12 that as we 
descend in energy the correlation length of a thermalized configuration 
becomes more and more dominated by the one of its inherent structure. 
We can guess that the divergence of the correlation length at the criti- 
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Variation of ^ CG with (5 



P 


Figure 4.11: The longitudinal and transversal Chiral Glass correlation lengths of the 
Inherent Structures in function of the inverse temperature j3 = 1/T. 
The trend is the same for both, specially when going towards lower 
temperatures and the derivative starts growing consistently. According 
to the most recent works /3 sg = 1/Tsg — 7.75. 


cal temperature is completely due to the underlying inherent structure 
scenario. 

In figure 4.13 we show the spin glass correlation function for the states 
thermalized at T = 0.19, and the same magnitude for their inherent 
structures. The correlation lengths are indeed very similar since the 
difference between them is indiscernible by eye. We can also appreciate 
the finite-size effects for the L — 64 samples, since it is quite clear that 
the correlation lengths do not go to zero. 
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4.4 Inherent structures at finite temperature 



Figure 4.12: The Spin Glass correlation length versus the inverse temperature /? = 
1/T. There are both the correlation lengths of the inherent structures 
reached with two different protocols, and of the thermal states. We see 
that although at infinite temperature (/? = 0) they are consistently differ¬ 
ent, when /? grows they become more and more similar, up to the point 
that one could argue that the divergence is due exclusively to the shape 
of the energy landscape. We see also how, lowering the temperature, the 
correlation lengths of the inherent structures depend less and less from 
the A used. 


Variety of the Inherent Structures and CG self-overlaps By 

comparing the energies obtained with SOR from T — oo with a very high 
A, with the ones obtained having first thermalized at finite temperature, 
but with a lower A, one could notice that it is possible to obtain the 
same energies in different ways. A naive thought could be that these 
are simply two different ways to achieve the same result. This assertion 
would be incorrect because the energy is not the only relevant observable, 
and as we already said in section 4.3, it does not characterize univocally 
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4.4 Inherent structures at finite temperature 



Figure 4.13: The correlation functions of the thermal state and its inherent structure 
at T = 0.19. The similarity between the correlation lengths is quite 
evident, as also the exponential trend. We see that at T = 0.19 the 
lattice size L = 64 is not enough since the correlation function does 
not go to zero (it is always far from being compatible with zero), while 
L = 128 is enough because the last ten distances are compatible with 
zero. 


a state. As a matter of fact magnetic susceptibilities, correlation lengths 
(fig. 4.14) and CG self-overlaps (fig. 4.16) may be very different for 
inherent structures wth the same energy, hence there is a whole variety 
of inherent structures for each energy. 

The CG self-overlaps deserve a separate discussion. It is commonly 
known that in Ising spin glasses there is no visual difference between a 
minimum and a random configuration, so that visual inspection can not 
discern one from the other. On the contrary, we noticed that this is not 
true in Heisenberg Spin Glasses. In fact, as one can see in figure 4.15, 







4.4 Inherent structures at finite temperature 



e < ls > 


Figure 4.14: The correlation lengths of the inherent structures, associated to their en¬ 
ergy. The dashed lines represent inherent structures from A = 100, 300 
protocol, starting from states thermalized at different temperatures. The 
solid line instead describes inherent structures from infinite temperature, 
with different A’s. If the energy were enough to characterize an Inher¬ 
ent Structure, each energy would have associated only one correlation 
length. On the contrary we see that same energies reached with different 
protocols have different £, so we can’t sample all the inherent structure 
population only performing SOR with different A’s from one tempera¬ 
ture, neither can we conforming ourselves to a single A from different 
temperatures. All the data comes from L = 128 lattices, except the 
points for A = 100,300 at T = 0.25,0.30,0.50, which can be recognized 
because of their lower correlation length, and the point for A = 10000 
and T = oo, whose study was so slow that we had to limit ourselves to 
few samples, which were obtained from L = 64 lattices. 


the CG self-overlaps are very different between thermal configurations 
and relative Inherent Structures. More specifically, we see that the CG 
self-overlaps are considerably smaller in the Inherent Structures, which 
have greater alignment properties. 
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Figure 4.15: The Chiral Glass self-overlaps of the inherent structures compared with 
the ones of the thermal states. The Inherent Structures have remarkable 
spin aligning properties, so, differently from Ising spin glasses, one could 
say whether a configuration can be a ground state or not, by simple 
visual inspection. 


Although it might seem, from figure 4.15, that the Inherent Struc¬ 
tures’ self-overlap is constant with temperature, one can correctly guess 
from figure 4.16 that those too descend appreciably with T. Conse¬ 
quently we can characterize an Inherent Structure also with its Chiral 
Glass self-overlap, and see, as it is well shown in that same figure, that 
this is another feature, other than the energy and the correlation length, 
that distinguishes different kinds of Inherent Structure. 
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4.5 Dynamics 



Figure 4.16: The CG self-overlap for Inherent Structures obtained with two different 
protocols. It appears clear that Inherent Structures with the same energy 
have different self-overlaps, therefore the energy does not individuate 
uirivocally the Inherent Structures. The error bars are greater in six of 
the points because they represent data coining from L = 64 lattices, 
differently from the others that come from L = 128 (we had to use 
L = 128 for lower T’s due to finite-size effects). Notice that lower T’s 
yield lower energies, so the same plot as a function of T would look 
similar. For this reason the reader can notice that the self-overlap of 
the inherent structure decreases with T, and is not constant as it may 
appear from figure 4.15 


4.5 Dynamics 

The self-correlation times (figure 4.17) yield more or less the same 
result obtained in ref. [8]. They diverge at a temperature significantly 
greater than T$c (or T C g), justifying the analogy with structural glasses. 
Also in this case we postpone the interpolations to the next section. In 
any case the divergence seems quite predictable also qualitatively, to a 
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T > Tsg — 0.12. Notice how the Spin Glass signal is sensibly stronger 
than the Chiral Glass. 



Figure 4.17: The two relaxation times Ts C aiar and tqg- The chiral sector carries 
information for much less time. 


4.6 Extrapolations 

Since most of the critical exponents of the 3-dimensional Heisenberg 
model are not known yet, we made some extrapolations that were possible 
with the data we collected. We have been able to estimate the critical 
exponents v and z, the exponent of the divergence predicted by Mode 
Coupling Theory (MCT), and we also wanted to give some estimates 
of Tqg by making extrapolations from far from Tsg, to show how large 
lattices far from the critical temperature are effective as small lattices 
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close to it, to make understand how it was possible to underestimate 
it so much to think it zero, and specially to remark that the nature 
of the overestimation of T c with the r’s is not a trivial misestimation 
effect due to being far from T S c, since that effect would result in an 
underestimation. 

4.6.1 An Estimate of Tsg 

We tried to estimate, the spin glass critical temperature by means of 
the divergence of the Inherent Structures’ correlation length. The scaling 
behavior near Tsg should be 

C 1 /v ~(T-Tsg)> (4.5) 

so if we plot versus (T — Tsg), and impose v to be the best 

estimate found until now, given in [17], i.e. v = 1.5, we can find Tsg from 
the linear fit’s intercept (figure 4.18). Except for the point at T = 0.5 
the data fits quite well, but yet we find T^T = 0.095(2), 20% lower 
than Tsg — 0.129 reported in [17]. This result was expected, since we 
are estimating Tsg with tools that are valid only near it, while our data 
are very far from it. Intuition leads to think that the data would assume 
curvature if we got closer to Tsg, until reaching a stable slope, so probably 
the real value of v greater than 1.5. 

Of course it is fairly reasonable that our estimate for Tsg didn’t repro¬ 
duce the results obtained in works who focused on the phase transition. 
This because the scaling laws are valid only close to Tsg , so for this aim it 
is worth to get some finite-size effects and work at smaller temperatures, 
around Tsg- 
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0 0.1 0.2 0.3 0.4 0.5 

T 

Figure 4.18: Power law fit of £, S q U imposing the value estimated in [17] v = 1.5. 

The point at highest temperature has been excluded from the fit since 
it didn’t respect the linear trend, and the scaling law is valid only close 
to Tsg■ The intersect gives us Tsg = 0.95(2), that underestimates it 
considerably. We see how using very large systems is not effective if we 
look for critical quantities. It is quite better to work closer to Tsg hr 
spite of the finite effects, since the scaling laws are valid only very close 

to T sg . 

4.6.2 An Estimate of the Exponent v 

From the same data set of figure 4.18 we can try to deduce the expo¬ 
nent v, by fixing Tsg = 0.129 as reported in [17], and fitting on v. We 
plot £ _1 versus T — T S c and expect a trend of the form 

r 1 OC (T - T SG y, (4.6) 

and obtain v = 1.146 ± 0.009 as it is shown in figure 4.19. If we 
compare it with the results in [17] we see that they get v{L = 8) = 
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4.6 Extrapolations 



Figure 4.19: Imposing the Tsg as the one evaluated in reference [17] ( Tsg = 0.129), 
we calculate the critical exponent v with a fit. The result v = 1.146(9) is 
compatible with the one obtained in [17] with smaller lattices, but closer 
to Tsg • They obtain v = 1.01 in a L = 8, and v = 1.35 with L = 12. 


1.01 ± 0.02 and u(L = 12) = 1.35 ± 0.05. Already L = 12 with measures 
taken close to Tsg gives a better estimate than ours. It appears clear how 
investigating the critical properties far from Tsg is equivalent to going 
very close to the phase transition with very small lattices, so are data is 
only good for an investigation of the pure paramagnetic phase. 

4.6.3 Self-correlation times 

This the most significant fit we present in this section, as it repre¬ 
sents one of the objectives of the present work. What we want to see is if 
the analogy with structural glasses holds and we can individuate a Mode 
Coupling transition, which predicts a power law diverging trend of the 
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relaxation times in the disordered phase, which is stopped, at lower tem¬ 
peratures, when we get in the region where activation leads the dynamics. 
The trend predicted by MCT is 

r oc (T - T c )~t. (4.7) 


100000 


100 t- 


10 


0.1 

T-T r fit 


^scalar r 1 
fit: A(T-T C P 


A = 0.453(8) 
C=3.33(2) 
T c =0.1489(8) 


Figure 4.20: Power law fit of the scaling of the scalar relaxation times Tg ca i ar in the 
paramagnetic phase. We encounter the same scaling, for 3-4 orders of 
magnitude, in supercooled liquids. The scaling individuates the diver¬ 
gence of the tS calar’s at T = 0.1489(8), which far in the paramagnetic 
phase, just as predicted by Mode Couping Theory for supercooled liq¬ 
uids. 


In figure 4.20 we can see how this prediction fits perfectly also for the 
EA Heisenberg model. We get a T c > T S g that overestimates the phase 
transition, just as MCT requests, and a perfect power law with exponent 
B = 3.33(2) which is of the same order of the those found by Kob in 
[36] (B = 2.5,2.6). We find T c = 0.1489(8), that is lower than the one 
found in [8] (T " c = 0.19(2)), but only two standard deviations away. We 
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4.6 Extrapolations 


stress that this T c shows itself as completely different from Tsg- In fact 
if it were different from Tsg just for the fact that we’re very far from 
it, our result would at least be consistent with the one of the previous 
paragraph, underestimating T S g , while we find T c » T S a » T <«. 

This result gives strong support to the possibility of a strong analogy 
between short-range interaction HSG and structural glasses. 

4.6.4 A Lower Bound for the Dynamic Exponent z 

Having data both for r and for £ we can think of trying to estimate 
the exponent z, that relates them through 

r(0 = if(l + 5r + ..-). (4.8) 

Unfortunately r(£) has strong scale corrections, since we only have 
£ < 8, so we only feel to give a lower bound for £. As we show in figure 
4.21 we can be reasonably sure that z > 4. but we cannot say how much 
it will grow. 
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Figure 4.21: The relation between relaxation time Tgcaiar and correlation length 
£Scalar follows a power law with exponent z close to the critical point. 
Since our data was far from the critical point it was not linear enough 
(in the logarithmic scale) to permit a decent fit, we limited ourselves 
to giving a lower bound for this exponent. The points connected by a 
solid line represent our measures, and the dashed line is the function 
/(£) = A£ 4 . It is clear how this function gives a loose lower bound of 
z>4. 


4.7 Minima or Saddles? 

To complete the frame that should indicate us if the analogy between 
the 3-dimensional EA Heisenberg model and structural glasses is too 
hasty we need to do a study of the density of negative directions in the 
hessian. In fact Mode Coupling Theory claims that the typical config¬ 
urations should be saddles until the temperature at which we have the 
topological transition. We need so to find if the typical configurations 
are closer to a minimum or to a saddle point. 
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4.7 Minima or Saddles? 


To this sake the algorithms we talked about until now are useless, 
since they are conceived to find minima of the potential energy and they 
would not make us perceive the presence of saddles. We need to find a 
new routine that returns us the stationary configurations mindless of if 
they are stable (minima) or unstable (saddles). 

The stationary configurations in our model are those for whom each 
spin sg is parallel to the local field hg = Yly^ 3-^=1 JxySy- The problem of 
looking for both minima and saddle points then is equivalent to minimiz¬ 
ing the spins’ component perpendicular to h, or maximizing the modulus 
of the parallel component. This way we could apply the algorithms we 
already used, to the energy functional 

H (saddle) _ _ J^(hg ■ Sgf (4.9) 

X 

and analyse the configurations out coming from this routine. 
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Conclusions 


We have made a systematic investigation of the behavior of the 3- 
dimensional EA Heisenberg model in the paramagnetic phase in the ther¬ 
modynamic limit, to see if we could remark features in common with 
structural glasses. Mode Coupling Theory predicts that the dynamics of 
structural glasses is dominated by the underlying structure of the energy 
landscape, so we made a deep study of the inherent structures of the 
model. 

For this sake the algorithms already used in literature were not enough 
effective, so we adapted the Successive Over Relaxation method to our 
problem and obtained a very performing routine. This algorithm de¬ 
pended on a parameter A, whose modification influenced highly its con¬ 
vergence time and the properties of the obtained inherent structures. We 
chose A who made our routine very fast, although it did not yield the 
lowest energies. This because we did not believe that the properties of 
the configurations we found were intrinsically different from the others, 
and also because this choice permitted us to maintain ourselves in the 
thermodynamic limit since lower energies were related to longer correla¬ 
tion lengths. We noticed also that the energy is not enough to identify 
univocally an inherent structure, since same energies found with differ¬ 
ent routines lead to configurations with completely different correlation 
lengths. 


81 



Conclusions 


We analysed the dependence of the energy of the inherent structures 
from temperature, and remarked a trend that was qualitatively similar, 
with almost no dependence at high temperatures, and a sudden decrease 
of their energies at lower temperature. Furthermore, we noticed that low¬ 
ering the temperature the correlation lengths of thermal configurations 
and inherent structures converge, so the divergence of the correlation 
length could be clue exclusively to the energy landscape. An analysis of 
the anharmonic term of the energy showed that its contribution is fairly 
small up to pretty high temperatures, so it is not so crazy to schematize 
the energy of a thermal state as harmonic agitation over a local minimum. 

We also noticed that the Inherent Structures of the Heisenberg Spin 
Glasses have particular alignment properties. In fact their Chiral Glass 
self-overlap is consistently lower than the one of the thermal states, and 
decreases with the energy of the Inherent Structure. This feature would 
be very helpful to recognize a Ground State, since it implies that, differ¬ 
ently from Ising Spin Glasses, the minima of the energy are qualitatively 
different from random configurations. 

Moreover, the data we had collected permitted us to make estimates 
for the critical exponents v and z of the model. Although they are not 
very reliable, since they come from states far from the critical temper¬ 
ature, where the scaling does not apply, they were useful to give some 
bounds that can be useful for further reference, since in literature there 
still are very few attempts of evaluating the critical exponents of the EA 
Heisenberg model. We focused on how working far from the phase tran¬ 
sition with very large lattices is equivalent to using small systems close 
to the phase transition, giving a measure of Tsg that underestimated it 
consistently, just as it happened when this study was approached using 
too small lattices. 

Later on we passed to the dynamic behavior of the model, since we 
wanted to see if we could observe the power law divergence of the rclax- 
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ation time in proximity of the dynamic transition, deep in the paramag¬ 
netic phase. Also this test had a positive feedback, since we observed the 
apparent divergence of the relaxation time for 4 decades just as in struc¬ 
tural glasses, the temperature at which the characteristic times would 
diverge was clearly greater than T$g and the scaling exponents of the 
same order of those obtained for supercooled liquids. 

All this data is consistent with the identification of the EA Heisenberg 
model as a possible lattice model with the same universality class of 
structural glasses. If this were true it would be very useful to study in 
an alternative way structural glasses. Further studies to make more solid 
this impression, would refer to the investigation of the configurational 
entropy of the model, and, work that is already in progress, a detailed 
analysis of the density of negative directions in the hessian matrix, to see 
if we can individuate the topological transition as its vanishing, with the 
arising of an activated dynamics. If those tests too will give a positive 
result it will be then enough evidence to start considering a rigorous, 
more engaged, theory of our conjectures. 
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Appendix 

Appendix 

A.l Fourier Correlation length 

We show how to get to the formula 3.37 for the correlation length, 
and show a computationally efficient way to calculate its terms. 

Given a generic held 0(F) 1 the action for a free Hamiltonian is 

S = ^ j d 3 x{(V(f)(x)) 2 + m 2 (j)(x)' 2 }, (A.l) 

from which descends the propagator 

<*<?> = S ■ * 2 =4l> 2 (t) (A.2) 

M=1 

where Z(fh 2 ) is an unknown function of the mass fh. We know that the 
propagator G is the Fourier Transform of the correlation function 

G(f) = + r) > ■ (A.3) 

X 

We can obtain the correlation length = 1/m 2 by making the ratio 
between the G calculated for two different impulses k. This way we 
would get rid of the factor Z{m 2 ). To this sake we use k — 0 and 

^or the purpose of our work <j> can represent both th Chiral Glass overlap and 
the Spin Glass overlap r a p. 
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k m in — ( 27 t/L, 0, 0), that is the minimum impulse with periodic boundary 
conditions, and we consider the ratio between 


F — G(k min ) 


Z(rh 2 ) 

4 sin 2 (-|) + fh 2 


and the magnetic susceptibility 


(A.4) 


X = G(0) 


Z(fh 2 ) 
fh 2 


(A.5) 


that yields directly 


c 2 = 1 [A _ p 

4 sin 2 (f) [ F J ' 


(A.6) 


Just to stay connected with the literature, we stress that with a finite dif¬ 
ferences approximation the correlation length in eq. A.6, can be written 
in differential form, since 


& = 


1 ,G(0) — G(k m i n ). 


k 2 

1X1 min 


G(k r 


dG 


G(k min ) dk 2 


\k=0 


(A.7) 


Calculating F and y An operative way of calculating F and y is 
made very simple thanks to the convolution theorem. Infact we see from 
equation A.3 that the correlation function is a convolution product of the 
field <fi with itself, hence its Fourier Transform G is equal to 

GW = k\m\ 2 , (A.8) 

where (j){k) is the Fourier Transform of the field cj)(x ). Remembering that 
we’re only interested in its value at k — 0, k m i n , we can use an explicitly 
computable form for 0( k) for F and y, helped by the fact that k m i n is 
parallel to the x-axis 

kkmin) = Y, ei2 ™ x/L Y ^(( n - ri y , n z )) = Y e^ /L P(f) (A.9) 

n x 'ft'yi'R'Z 'kf'X 

where n x , n y and n z are the components of the position vector x, and we 
recognize the plane variables P(x) defined in section 3.3.3. 



A.2 Integral Estimators 


We now have a numerically simple way to express F, from which the 
one for x descends trivially, 

x = i(X>W) 2 < A - 10 ) 

X 

F = (A.11) 

Tlx Hy i'Fjz 


A.2 Integral Estimators 


We want to show that in the thermodynamic limit the integral esti¬ 
mators ^ 02 - ) are exactly equal to the correlation length We remind 
its definition 


o (02) = 



■h 


L/2 


E ' , 2 C(r) 


(A.12) 


The minimum impulse propagator, keeping in mind that we work 
with periodic boundary conditions, can be rewritten as 



r =1 


(A.13) 


where the sum on the sines disappears because we’re summing an odd 
function over a symmetric interval. If we pass to the large L limit, in 
which we want to demonstrate the identity, we can develop the cosine 
for small angles and neglect the term C (^), since C(r ) vanishes at long 
distance. Hence 


L/2—1 


F = d(0)+2^C(r)[l- —+ o(-)| = 

r=l 


L 2 


r 4 

D 


= x 


A) 2 L + o(^). 


L 4 ' 


We can insert this expression in the formula for getting 

1 r 1 




W(f) L l-(f) 2 J + o(^) 


1 ], 


(A. 14) 
(A.15) 


(A.16) 
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from which if we develop another time for large L, keeping in mind that 
for small x we have ~ 1 + x, and take the limit L —> oo we obtain 

1—X 

finally 

lim £f= ,/^ = £ (02) . (A.17) 

L ^°° y X 

Figure A.l shows how the correlation lengths and that differ 

quite much for small system size, converge as the system approaches the 
thermodynamic limit. 


Spin Glass Correlation Length 



1/L 


Figure A.l: Correlation lengths of configurations obtained with 3 different protocols 
(A = 10,100,1000). We see how in small lattices, far from the thermo¬ 
dynamic limit, and belonging to the same type of configuration 
are very different, since their equality is valid only in the thermodynamic 
limit. With growing system sizes the two correlation lengths start con¬ 
verging considerably. The convergence of these two correlation functions 
can be used as index to establish if the thermodynamic limit is valid or 
not. 
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A.3 Asymptotic values of Binder’s Cumu¬ 
lants 


We show how to extract the asymptotic value B^ of the SG Binder 
Cumulant defined in section 3.3.5 in the case of configurations at 
infinite temperature (random configurations). The calculation of B is 
analogous so it won’t be shown. We remind the definition 

SG _ < tr(QQ^) > 2 = < q 2 > 2 
L < tr(QQi ) 2 > <q A > 

and proceed calculating separatly numerator and denominator for spins 
with m components. We start with the numerator 


< tr(QQ ] ) > 


^ ^ ' Qa 

a,(3 


* >= 


< 44444444 >= 


^ 1 ^ 1 0 x,a°x,l3 0 y,aL 0 y,j3 

a,0 x,y 


E 

x,y 

E« 


< (4“> - 4 a) ) >< (4 6 > • 4 6 >) >= 


x, y 

N 


and proceed with the denominator, where due to the square we have two 
more position indexes 


< tr(QQ ] f > = y <(4 a) -4 a) )(4 6) -4 6) )(4 a) -4 a) )(4 6) -4 6} ) 

x,y,z,t 

= y <& )4 € 4> + 

x 

+ Y < 4 a)2 4 ft)2 4 a)2 4 6)2 > + 

x,z 

+ ?Y< (4 0) • 4 a) ) 2 >< (4° • 4 b) f > • 


>= 


x,y 


We compute the average < (s x • s y ) 2 > by putting ourselves in the refer¬ 
ence frame of s x , aligning to it our first axis. This way the scalar product 



A.3 Asymptotic values of Binder’s Cumulants 
is equal to the first component of s y . Hence 

< (fix • Sy) >= Sy,l = — 

since X1 q=i ^ i = 1 and there is no preferential direction. Therefore, 
summing up all terms we have 

< tr(QQ t ) 2 >= N + N(N - 1) + \N(N - 1). (A.18) 

rrr 

Recomposing the cumulant and performing the thermodynamic limit we 
get 

B f i + 2 ( A . 19) 

m 

So in our case m — 3 we have = y = 1.2. 
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